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ABSTRACT 


This thesis devoted to studying the properties of a digital M-ary signaling scheme called 
MFSK/QPSK which combines the features of M-ary FSK and QPSK modulation. This 
scheme is shown to be a form of biorthogonal signaling, and a receiver that is optimum 
in error rate sense is introduced and analyzed in order to obtain its bit error rate (BER) 
performance when signal interference can be modeled as additive white Gaussian noise. 
The bandwidth of the MFSK/QPSK signaling scheme is derived and compared with that 
of M-ary PSK and M-ary FSK signaling. MFSK/QPSK signaling results in receivers 
that provide error rates much lower than corresponding receivers for M-ary PSK and 
slightly lower than those for M-ary FSK, while maintaining a bandwidth efficiency at 
least three times greater than that of M-ary FSK though not as great as that of M-ary 
PSK signaling. A significant effort has been devoted to examining a detection method 
known as direct bit detection (DBD) which has been applied to M-ary PSK and 
MFSK’QPSK signaling. DBD receivers are introduced, and the BER of the 
MFSK/QPSK DBD receiver is derived and compared with that of the optimum (symbol 
discrimination) receiver. Lastly, the effects on performance due to Rayleigh fading are 
analvzed and evaluated for M-ary PSK DBD receivers as well as for the MFSK’QPSK 


optimum (symbol discrimination) and DBD receiver. 


Led NERODUCTION ..... 6500.00 0 st ee gees - eNSER,.. Sanne ] 
A. .BACKGROUND 2.000... SER emt ] 
B. CONTENTS AND SUMMARY OF RES U EiiSa = ee 1 
C. M-ARY PHASE-SHIFT KEY IiiGee ee ee 2 
D. M-ARY FREQUENCY -SHIP IKE aN G Sr eee 3 
E. RECEIVERS*OFsMcARY PSK ANDEM-ARY FSke... a... 5) eee 3 
F. SYMBOL ERROR RATES FOR M-ARY PSK AND M-ARY FSK ...... 7 
G. BANDWIDTH REQUIREMENTS ©. ..oeercne 10 

Il. DIRECT BIT DETECTION METHODS FOR "I-AnyY PSK ee 12 
A. "M-AIRY PSK DIRECT BIT DETEGRION KECEM ERS... oe 2 
B. PERFORMANCE OF DBDSMECEIVERS FOR M-Atwy PSK 3 eee }2 

TIT. MPFSK?OPSK yoo. o.oo 20 
A. MFSK/QPSK AS A BIORTHOGONAL SIGNAL SE) eee 20 
B. OPTIMUM RECEIVER .... 0.23 ee 27 
C. OPTIMUM RECEIVER BE Rae eee ee 38 
D. DBD RECEIVER FOR, MIESK OES 0) 2 4] 
E. MFSK;QPSK DBD RECEIVE R@O 2 a ee eeniene . cgeerenee 43 
F. BANDWIDTH OF MESK (ORS igen gre ere 46 

IV. RAYLEIGH PADDING (2. See ete eee rere ne sere 50 
A. THE RAYLEIGH FADING MOD ee. ee 50 
B. EFFECT ON BER OF MARY PSK USING DED RECEIVER. ...2aeee 50 
Cc. EFFECT ON BER OF MFSK/ORSK USING DED KE GEIVER |e 61 
D. EFFECT ON BER OF MESK/QPSK USING THE OP TIME MicE- 


CEIVER on che gee 0 6 wg cee ne ee a 62 
V. RESULTS 6.06 ca co. +s RURRRRRRR en epee iretee tat rane er 64 
A. PERFORMANCE COMPARISONS "322... 5 505 64 
B. BANDWIDTH EFFICIENCY 22253 ee 70 


ee sore eG FADING «2.6.6. ee es 70 


Sl CUO NICILT OS IOO S oh ne 6566. 7 
Be SIDIOR ay NS Te OR US a 78 
eS Oe ees OPSK BER COMPUTER PROGRAM ............ 82 


APPENDIX  C. 16-PSK RAYLEIGH FADING BER COMPUTER PROGRAM 85 


APPENDIX OD. MFSK/QPSK RAYLEIGH FADING BER COMPUTER 


Oo CPE 9 gg os ooo 6 Ca 00 2 a a 
CLS IE RIE RE IIE SSC) SR er a gB) 
SLE AIL ITER TSEC ON S130 96 


Table 
Table 


Table 
Table 


Table 


LIST ‘OF TABS 


l. FUNCTION DEPINIVIONS 2322353 ee 18 
2. PERFORMANCE COMPARISONS AT KEY SYMBOL ERROR 
RATES oc os s:ishig: aie micelle adage oe cpm eee 64 


3. PERFORMANCE COMPARISONS OF MESK/OQPSK RECEIVERS as 
4. BANDWIDTH-EFFICIENCIES FOR M-ARY PSK, M-ARY FSkK; 

AND MPSKIOPSE “0s seas cee eS eens og) er 70 
5. EFFECT OF RAYLEIGH FADING USING VARIOUS RECEIVERS 


wi 


Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 


Figure 


Oo COO HN mH BW YN — 


10. 


LIST OF FIGURES 


So emem@onstellanons fora) QPSK and b)8-PSK ...............- 
2 IZ OGNTIOS S CECIONIG CIO) ES) GIR a 
Me dem @eropecinimmamveaaty FSK 2... ee ee ee 
MiitcerivemomuctiretOmM-ary PSK... 1. ds i. see ee ee es 
Mel eecmveIm CMTC NOI vi-aly ESK 2.64. aca ee ee ee ee ee 
mplectemiconaiy othe s-Fok DBD Receiver .........60s000055- 
meee miapramon tie 16-PoK IBD Receiver ........2...005505. 
. Signal Constellation for 8-PSK Showing Variable Phase Parameters 

. Signal Constellation for 16-PSK Showing Variable Phase Parameters . 
P eanomincmnesemiamoneon JS (OPSK 9... 6 a ee 
Pe DiacmamlrOlcanencialCOnmelauion KECEIVED «22. ee eee ee eee 
me Picram on @orrelation Receiver for 2FSK/QPSK .... . = sags aes: - 
me Pincwinsonm bh Oekecciver hor 2FSK/OPSK 2.............26...4. 
Me camudcespectnum for 2ESK/OPSK « . wnuiaa.s same s sc eew ss oe 
me iacmecemspecthunitor MESK-OPSK 2... cee ee ee ee 
mE Siwe@rsic Bit Error Performance {k = 3 bits’symbol) .......... 
me ok OPSK Bit Error Performance (k = 4 bits/svymbol) .......... 
moro Poke Pi Emor Penornmance (k — 5S bits’simbol) .......... 
melon sheOrsk bit Error Performance (k — 6 bits symbol) ......... 
. 8-PSK DBD Receiver Bit Error Performance in Rayleigh Fading .... 
. 16-PSK DBD Receiver Bit Error Performance in Rayleigh Fading 
22. 
23: 


MFSK/QPSK DBD Revr. Bit Error Performance in Rayleigh Fading 
MFSK/QPSK Opt. Revr. Bit Error Performance in Rayleigh Fading 


Vil 


— ~) ~~!) 
pa 


om 


~) 
Gd bo 








I. INTRODUCTION 


A. BACKGROUND 

There exists an increasing need for high-data-rate, reliable digital communication 
systems. The reasons for this need stem from the desire by an increased number of users 
to transmit larger amounts of data via digital communication schemes which provide 
more reliability, flexibility, and lower cost when compared with analog communication 
systems. These increased demands come at a time when the available bandwidth in the 
free space channel is becoming a limited resource. These data transmission requirements 
coupled with limited channel resources are the driving force behind the continual search 
for digital modulation techniques which operate at low error rates and are bandwidth 
efficient. Present signaling schemes in use are the Well known M-ary Phase Shift Keying 
(M-PSK) and M-ary Frequency Shift Keving (M-FSK). These signaling schemes are 
briefly described in Sections C and D of this chapter, and the corresponding receiver 
structures for recovering the digital information are described in Section E, while the 
performance of these receivers, specified in terms of the symbol error rates, for both 
M-ary PSK and M-ary FSK modulation are given in Section F. The bandwidth re- 


quirements associated with these two modulation schemes are given in Section G. 


B. CONTENTS AND SUMMARY OF RESULTS 

This thesis explores a hybrid form of digital modulation, namely MFSK/QPSK, 
which embodies features of both M-ary FSK and Quaternary (or four phase) PSK. The 
analysis compares and contrasts the properties of this modulation scheme to the widely 
known M-ary PSK and M-ary FSK modulation. In Chapter 3, MFSK, QPSK is shown 
to be a form of biorthogonal signaling, and therefore its corresponding receiver per- 
formance is easily revealed given what is already known for this class of signals. The 
optimum MFSK/QPSK receiver that recovers the digital information with minimum 
probability of error under specific noise interference assumptions is analyzed from a 
performance standpoint, which traditionally has been given in terms of the Symbol Error 
Rate (SER). Because a more meaningful performance parameter is the receiver's Bit 
Error Rate (BER), Chapter 3 has been devoted to deriving the receiver's BER as a 
function of signal-to-noise ratio, along with the bandwidth requirements and the band- 
width efficiency of the MFSK/QPSK scheme. Comparisons are made between the 


MFSK, QPSK scheme and conventional phase and frequency modulation methods with 


respect to SER performance and bandwidth efficiency in Chapter 5. The results show 
that the MFSK/QPSK scheme when compared to M-arv PSK, does not attain the high 
bandwidth efficiencies achieved by M-ary PSK, yet outperforms M-ary PSK substan- 
tially in terms of SER performance, while performing slightly better than M-ary FSK in 
terms of SER performance and only requiring one third to one fourth the bandwidth 
required by an M-ary FSK system. 

Additionally, this thesis addresses the properties of a class of digital communication 
receivers Which decode the individual bits that make up a symbol directly, without need 
for symbol-to-bit mapping post-detection circuitry. These are referred to as direct bit 
detection (DBD) receivers, because each of the bits that make up a symbol are deter- 
mined independently. These receivers have been proposed and analyzed for M-ary PSK 
Ref. | [pp. 461-463] as well as for the MFSK/QPSK signaling scheme Ref. 2 [pp. 16-101] 
under discussion here. DBD receiver structures for M-ary PSK are introduced in 
Chapter 2, and previously derived expressions for their BER performance are presented. 
In Chapter 3, a DBD receiver structure for MFSK/QPSK signaling is described, and 
analysis is carried out in order to derive its BER performance, with corresponding results 
presented in graphical form in Chapter 5. Also in Chapter 5 the BER performance of 
the MFSK,/QPSK DBD receiver is compared to that of the optimum receiver that per- 
forms symbol detection followed by symbol-to-bit mappings prior to delivering binary 
data to the intended user. Furthermore, the results demonstrate that while the DBD 
receiver does not perform as well as the optimum receiver in BER sense, its simplified 
receiver structure and direct bit detection properties offer other potential advantages. 

In choosing a modulation method, it is important that the effect of svstem disturb- 
ances such as Ravleigh fading, which is common in digital microwave radio transmission 
applications, be analyzed in order to fully understand all aspects of expected system 
performance for the different signaling schemes under consideration. The effect of 
Rayleigh fading on the BER performance of the MFSK/QPSK receivers as well as of the 
M-ary PSK DBD receivers is analyzed and evaluated in Chapter 4. The results demon- 


Strate that in all cases, Rayleigh fading can severely degrade receiver performance. 


C. M-ARY PHASE-SHIFT KEYING 

M-ary signaling techniques involve forming groups of k bits to generate M = 2 
symbols which are assigned to M different waveforms suitable for transmission over a 
channel. One form of M-ary signaling, namely M-ary Phase-Shift Keving (M-PSK), is 


a bandwidth efficient digital modulation scheme. The phase of the transmitted 


sinusoidal carrier waveform can take on one of M different phases to represent M dif- 


ferent symbols. The M-ary PSK waveform can be described mathematically as 


s{1) = je cos| 2afge + AEE | (1.1) 


fori = 0, l,..., M-1 and 0<7<T,, where generally (as is the case here) the phase dif- 
ference between waveform pairs is taken to be constant. Each symbol or waveform 
${t) is transmitted at a rate of R, = al symbols/sec, where 7, is the time duration of each 
waveform s{t). For k=2 and M=4, the resulting modulation form is known as 
Quaternary Phase-Shift Keying (QPSK), where the phase difference between waveform 
pairs is radians. For 8-PSK, in which k = 3 and M = 8, waveform pairs are separated 
by a phase oie radians. Signal constellations for QPSK and 8-PSK as well as symbol 
assignments for each signal (or phasor representation) are presented in Figure | on page 
4. 


D. M-ARY FREQUENCY-SHIFT KEYING 
M-arv Frequency-Shift Keying (M-FSK) is another form of M-arv signaling, in 


which the waveforms are mathematically expressed as 


r= joe cos 2nft (1.2) 


where the frequencies f, typically are given by 
Si=So + iAf Ch.3) 


fori = 0, 1,.., M-1 andQ<?r<T,. The signal duration is normally chosen so that fT, 
is an integer or at least f7,>> 1 holds. The frequency separation Af is often set equal 
to R, providing minimum frequency spacing between signals s(z) , while forming a set 
of orthogonal signals. For k=2 and M=4, the resulting modulation is Quadrature 
Frequency-Shift Keying (QFSK). Magnitude spectra for QFSK and M-ary FSK ( for 


arbitrary M ) are shown in Figure 2 on page 5 and in Figure 3 on page 6, respectively. 


E. RECEIVERS OF M-ARY PSK AND M-ARY FSK 

Signals transmitted through a physical channel always encounter interference. The 
SER results given later in this chapter and the BER analysis presented in Chapter 3 are 
obtained under the assumption that this interference can be modeled as a sample func- 


tion of an additive white Gaussian noise (AWGN) process having double sided power 


sin(2 mf;t) 


01 


00 
a) cos(2 fot) 


sin(2 nf,t) 


b) cos(27f ot) 





Figure 1. Signal Constellations for a) QPSK and b) 8-PSK 
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spectral densitv (PSD) of 2 Under this assumption, receivers can be derived that are 
optimum in the sense that the probability of a symbol error is minimized. One structural 
form of an optimum receiver for M-ary PSK is shown in Figure 4 on page 8. This re- 
ceiver as Well as all receivers considered herein is assumed to be coherent, meaning that 
the phase of the sinusoidal carrier is known at the receiver. The concept of coherence 
applies to this receiver in that the cosine and sine functions which are correlated with the 
received signal are in phase with the unmodulated carrier sinusoid. Because the receiver 
correlates the received signal with a set of sinusoidal functions, which in this case con- 
sists of each of the possible transmitted signals, it is referred to as a correlation receiver. 

An optimum receiver structure for M-ary FSK is shown in Figure 5 on page 9, 
which also corresponds to a correlation receiver. The correlator with a local signal input 
which most closely matches the incoming signal, will produce in the absence of interfer- 
ence, the largest output so that the receiver can decide which signal was sent according 


to which correlator output is the largest. 


F. SYMBOL ERROR RATES FOR M-ARY PSK AND M-ARY FSK 
For coherent demodulation of M-ary PSK using the optimum receiver shown in 
Figure 4, the SER is given by Ref. 3 [pp. 204-207] 


ntl ee oe PIS ees 
[as a aay, +0( Aaa: 





ez (1.4) 
— Zz [x * exp( —}) | — 20(y2 cot ) Jay 
va 70 
where 
o)= [Pea (1.5) 
g 


The integral on the right hand side of Eq. 1.4 has no known closed form, and therefore 
must be evaluated numerically if plotting of P,{e} as a function of = is desired. For 
INg 
high values of, the SER can be tightly upperbounded by 
4X 


Pee city 
Pfey< 20 y; sin 47 (1.6) 
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Figure 4. _- Receiver Structure for M-ary PSK 
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Figure 5. Receiver Structure for M-ary FSK 


For M-ary FSK modulation and an AWGN interference model, the optimum co- 
herent receiver of Figure 5 on page 9 achieves a probability of correctly decoding a 
symbol given by Ref. 3 {pp. 200-203] 


sad ae M-1 
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where o = > and g(f) is the zero mean unit variance Gaussian distribution, or 





This expression can be rewritten as 


rio-| [!-o(--VE)]sow-1-ra a 


Also, the SER can be closely upperbounded for high values of = by 
0 


| TE, 
P(e} < (M — no(./ a (1.10) 


G. BANDWIDTH REQUIREMENTS 
For M-ary PSK, the carrier is switched at the rate of R, times per second. Since 


there are k bits associated with each symbol, the bit rate R, can be expressed as 
R,=kR, (Itsy 
The null-to-null bandwidth W for M-ary PSK is given by 


2R 
W = 2R, => (1.12) 
For M-ary FSK the carrier is also switched at the symbol rate R,. Looking at Fig- 
ure 3 it can be seen that the total null-to-null bandwidth left-most main lobe to right- 


most main lobe for M-ary FSK can be written as 
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(1.13) 
From Eq. 1.13 it is apparent that the null-to-null bandwidth of M-ary FSK increases 
with increasing values of k, however Eq. 1.12 shows that the null-to-null bandwidth of 
M-ary PSK decreases with increasing k. This property makes M-ary PSK a desirable 
modulation method when bandwidth must be minimized. In the next chapter M-ary 
PSK receivers of a different structure than the one shown in Figure 4 are introduced, 


and their bit error rates presented. 


Il. DIRECT BIT DETECTION METHODS FOR M-ARY PSK 


A. M-ARY PSK DIRECT BIT DETECTION RECEIVERS 

M-arv PSK modulation as discussed in Chapter 1, Section A, is a signaling scheme 
that achieves bandwidth conservation at the expense of signal-to-noise ratio required to 
maintain a certain bit error rate performance level. The two forms of M-ary PSK sig- 
naling discussed in this chapter are 8-PSK and 16-PSK. For 8-PSK, k= 3 bits/symbol 
and 'f=8 symbols, while for 16-PSK, k =4 bits/symbol and 4f=16 symbols. The 
M-ary PSK receiver discussed in Chapter 1, Section E and most other optimum receiver 
structures discussed in the literature, require symbol-to-bit mapping to be carried out 
after a decision has been made as to what signal was transmitted. A receiver structure 
which does not require this svmbol-to-bit mapping is proposed in Ref. 4 [pp. 3-5] and in 
Ref. 1 [pp. 461-463] for 8-PSK and 16-PSK, respectively. Block diagrams for these re- 
ceivers are shown in Figure 6 on page 13 for 8-PSK, and in Figure 7 on page 14 for 
16-PSK. Because each of the bits are determined independently by these receivers by 
appropriate processing of the correlator outputs, these receivers are referred to as direct 
bit detection (DBD) receivers. The processes used by DBD receivers for recovering bits 
make it possible to adjust the BER of the individual bits within a symbol by means of 
varving the angles between the signal waveforms. These angle adjustments are shown 
in Figure 8 on page I5 for 8-PSK, and in Figure 9 on page 16 for 16-PSK, where o and 


B are the adjustable phase parameters. 


B. PERFORMANCE OF DBD RECEIVERS FOR M-ARY PSK 

A BER analysis for the 8-PSK DBD receiver was carried out by Ref. 4 [pp. 5-14] in 
which a coherent receiver was assumed, and AWGN used as the interference model. 
The three bits associated with each symbol are identified as 1B, 2B and 3B corresponding 
to the most significant bit, the second-most significant bit and the least significant bit, 
respectively. The results of the BER analysis showed that the individual bit error prob- 


abilities were not equal for the 1B, 2B and 3B. In fact these probabilities are given by 
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Figure 8. Signal Constellation for 8-PSK Showing Variable Phase Parameters 


sin(2 rf) 


s = 0011 s,= 1011 
SF 0001 5S 1001 
Q 
, = 0101 S = 110m 
s = 0111 s =n 
ZB . 
Ya") cos(27fot) 
s, = 0110 $= 1110 
5 = 0100 d Ss 1100 
s = 0000 = 1000 
10 = 0010 s= 1010 


Figure 9. Signal Constellation for 16-PSK Showing Variable Phase Parameters 
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The overall BER for 8-PSK using the DBD receiver can be expressed as 
Pie = = Pr{1B incorrect} +> Pr{2B incorrect} (2.3) 


For 16-PSK, there are four bits associated with each symbol which are identified as 
1B, 2B, 3B and 4B corresponding to the most significant bit, the second-most significant 
bit, the third-most significant bit, and the least significant bit, respectively. The error 
probability associated with each of these bits is different as given by Ref. 1 [pp. 461-463], 
under the assumption of a coherent receiver and an AWGN interference model. The 


bit error probabilities for the individual bits in 16-PSK are given by 
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where the integral in Eq. 2.6 has not been solved analvtically, but has been evaluated 
numerically in Ref. 1 [pp. 466-467], and the functions f([, «, B) and g{I, «, B) are de- 
fined in Table 1. 





Table 1. FUNCTION DEFINITIONS 
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The overall BER for 16-PSK can then be expressed as 


P,fe} = 


> Pr{1B incorrect} + = Pr{2B incorrect} + > Pr{3B incorrect} 


(2.7) 


These results have been presented here so that results in the sequel can be compared to 


the BER results associated with DBD coherent receivers for 8-PSK and 16-PSK signal- 
Ing. 


Ill. MFSK/QPSK 


A. MFSK/QPSK AS A BIORTHOGONAL SIGNAL SET 

If the features of M-ary FSK and QPSK modulation are combined, a hybrid digital 
modulation scheme results where signals are both phase modulated and frequency 
shifted. This new modulation technique will have four signals at each of the operating 
frequencies corresponding to the four different possible phase modulations. The signal 


set can be expressed as 








s(t) = ae cos2nfir + iS] (3.1) 
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fori = 0) 1, 2,3 and0 <ie7, , and 
50) = 2a eee (3.2) 
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fori = 4, 5,6, 7 andQ 271-7 


Integrating the square of the signal s,(r) over the interval 0 <7 < T, yields 
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If the product /,/, is an integer, or iff,7 = = | einen 
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which corresponds to the signal energy. Similar integration yields identical results for 
S(E) alee, aan 
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Applying the Gram-Schmidt orthonormalization procedure Ref. 5 [pp. 266-269] to 
the signal set described by Eqs. 3.1] and 3.2 generates a set of orthonormal basis func- 


tions {@,(7)} having the property 


T, m 1; for i=Jj 
i onodjoa =| ay i (3.5) 


so that the signals sr) can be synthesized using a linear combination of basis functions, 
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s(t) =) sybft)i §=0,1, ...,7 (3.6) 
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where the s, coefficients are defined as 


T, 
= 5(1)@ (ode (Gap) 


The first basis function is obtained by normalizing s,(t). From Eq. 3.4, we have 


< es) = cos[2nfit] (3.8) 


so that s,(2) = ./E ¢,(0), and the coefficient s,, = /E . From Eq. 3.7, 5, is given bv 
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Defining the auxiliary function y,(r) as 


WD = 51(2) — 54,6119 3 


= 5,(7) 


Where the second equality follows from Eq. 3.9, the second basis function can be ob- 


tained by normalizing the auxiliary function. That is 


$= =e le cos(2nft +>) (3.11) 


Using next the coefficients s,, and s,, to obtain the auxiliary function y,(r) as 
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results in Eq. 3.12 becoming equal to 


v(t) = (0) + VE $,(9) 


= lee cos(2afittn)+_/ = cos(2nft) (3.15) 


e ae [ — cos(2af;t) + cos(2a/,2)] = 0 





This indicates that s,(t) is completely described as 


s() =—- VE $,(¢) 


and thus no new basis function is introduced by s,{1) . 


The coefficients for s,(t) are given by 
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so that the auxiliarv function y,(z) becomes 


val ) = 83(2) — 532(2) 
= = cos(2mfjr + 3) + / oe cos(2nfyt +>) 


_ = [ — cos(2afit ++) + cos(2afr +4 )] =0 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


Thus, as with s,(z), s,(t) introduces no new basis functions. Continuing this procedure 


for the signal s,(r) , we have 


23 


sy] sald ae 


T, 
-| am cos(2zf,t) in cos(2nf\1)dt 
0 


25 


Ege 
= | cos(2zf52) cos(2zf 1dr = 0 
Ts Jo 





(3.20) 


where the last equality is due to the assumption that the frequencies f, and f, are spread 


far enough apart to maintain orthogonality as in binary FSK modulation. Similarly 
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so that the auxiliary function y,(z) becomes 
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yalt) = salt) — ) Sy 00) = 5360) 


J=1 


Normalizing y,(z) in order to obtain the third basis function, results in 


$3() = - —_ ee eos an) 


Next, evaluating s, for } = I, 2,3 yields 
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and 
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0 


2 
=— | cos(2zf,t + Zz) cos(2zf,t + = )at =0 
i” ps 2 





where again, the frequencies f, and f, are chosen to maintain orthogonality as explained 


above. Furthermore, 
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so that with all three of the coefficients equal to Zero, s.(r) must be orthogonal to all of 
the basis functions generated thus far. The fourth basis function is then obtained by 


normalizing s.(z) , that is 





g(t) = “~ = Slee cos(2nfzt + >) (G27) 


It can be easily shown that the remaining signals s,(r) and s.(r) can be expressed as 


s(t) = — VE $3(0) 


i 3.28 
s(t) = — JE $4(2) sa 


Thus the coefficients s,, and s., are both equal to — JE. Since the set {,(t)} represents 
an orthonormal basis in a vector space, the signals s() can now be interpreted geomet- 


rically as vectors having components equal to the s,, namely 


[Oe 5 =, Sere =0 te 1 (3.29) 


where N is the number of linearly independent orthormal basis vectors, and M is the 


number of signal vectors.! It is well Known that 


if 
| s;(t)dt = > =|5\?=6 (3.30) 
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where E, is the energy of the ith signal. For the case of 2FSK/QPSK, M is equal to 8, 


N is equal to 4, and the vectors representing the signals are 


% =(+VJE, 0, 0, 0) 
3,=(0, +JE, 0, 0) 
S=(— JE e000) 
s,=(0, — VE, 0, 0) 
3,=(0, 0, +/E, 0) 
0D. 0 ae) 
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5 = (0090, — 9 5n) 
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(3.31) 


Since there are eight different possible transmitted signals, each signal can represent 
three bits of information. 
The average energy and also the dimensionality of the signal set can be minimized, 


by subtracting from each signal a constant vector a such that 


M-1 
> PS 5 — al? (3.32) 
i=0 
is minimized, where P{s} is the a priori probability than the signal s(t) has been trans- 
mitted. Such vector a can be obtained from 
M-1 ; 
a=) PS}5= EG) Gea) 


f=) 


1 The general practice of labeling the number of signals as M is followed here, but the letter 
M in “MFSK;QPSK” represents the number of frequency shifts. 
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For a signal set of the type described above, assuming all a priori probabilities P{s,} are 


equal, Eq. 3.33 for a reduces to 


M1 
& =0 (3.34) 
i=0 


so that there is no minimum energy transformation that would reduce the dimensionality 
or the average energy of the signal set. 

Figure 10 on page 28 shows the vector representation of the 2FSK/QPSK signal 
set shown as eight vectors in a four dimensional vector space. A signal set consisting 
of N orthogonal signals augmented with the negative of each signal is referred to as a 
biorthogonal signal set. As shown by Figure 10, the 2FSK/QPSK signals form this type 
of set. For a general biorthogonal signal set, the number of dimensions of the vector 


space is always equal to half the number of signal vectors. That 1s 
M=2N (3539) 
while the number k of bits represented by each signal is given by 
k=log,M (3°36) 


Finally, the number K, of signal frequencies in the MFSK/QPSK scheme is one-fourth 


the number of signals or 


M=4K, (3.37) 


The best assignment of bits to the M signals of the biorthogonal signal set is to 
match those signals separated by the greatest Euclidean distance in the vector space, 
with symbols having the greatest Hamming distance. With this rule in mind the bit as- 
signments become obvious. That is, signals which are antipodal are assigned symbols 
which correspond to binary complements. Once this assignment is done, all other vec- 
tors’ Euclidean distances are equal, so that any other assignment patterns can be freely 


chosen. A typical bit assignment is shown for the case of M = 8 in Figure 10. 


B. OPTIMUM RECEIVER 
With the basis functions having been specified, available analvtical results on the 


structure and performance of receivers for discriminating biorthogonal signals received 
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Figure 10. Vector Representation of 2FSK/QPSK 
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in additive white Gaussian noise make it possible to design an optimum receiver for the 
discrimination of MFSK/QPSK signals. The received signal r(t) is assumed to consist 
of a signal component s,(1) plus a sample function of a white Gaussian noise process n(t). 
That is, 


r(t) = s{t) + n(0) (3.38) 


As shown by Ref. 5 [pp. 211-235] the optimum receiver takes on the structure of a cor- 
relation receiver. A block diagram of the general correlation receiver is shown in 
Figure 1] on page 30. The signal r(t) is processed by a bank of correlators, where the 
number of correlators in the receiver is equal to N, the number of dimensions in the 
signal vector space model. The effect of the additive white Gaussian noise is to add to 
the signal vector a noise vector n having N on components, each of which is 
a zero mean, Gaussian random variable with variance —— . Thus the correlator outputs 


2 
r, form the components of a random vector r, i.€:, 


T=s5+n=(r, m, .--, Iy) (3.39) 


After generating the vector r , the components r, are processed through a weighting 


matrix to generate the dot products d, which are given by 


=Te 5 = pie (3.40) 


j=l 


The d, are then summed with bias terms c, where 


=> (No Inf Pom} - ele r= oe veel (3.41) 


Finally, all summer outputs I, are routed to circuitry which decides which signal was sent 


according to the rule: 


jj i ee VY k=@) l, = el cag 


: : ; (3.42) 
Then decide st) was the transmitted signal. 


For the signal vectors of Eq. 3.31, assuming P{s(r)} = a for all /, the structure of 
the optimum receiver simplifies from that of the general correlation receiver to the 
2FSK/QPSK receiver shown in Figure 12 on page 32. The 2FSK/QPSK scheme can be 


2g 





Figure 11. | Diagram of a General Correlation Receiver 
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extended bevond the case considered by increasing the number of signal frequencies. 
For example, in order to increase the number of bits per signal from three to four, the 
number of signals M must be doubled to 16. This is accomplished by doubling the 
number of frequencies K, from two to four, and thus doubling the number of basis 
functions from four to eight. This change will be reflected in the receiver structure by 
a doubling of the number of correlators to eight. The new signal set is also 
biorthogonal, and further extensions to MFSK/QPSK result in N basis functions, K; 
frequencies, M signals, and k bits per signal. The formula for the probability of a correct 


symbol decision for the biorthogonal set is given by Ref. 5 [pp. 261-263] as 


N-1 





P{C} = | ila — JE vs —2 | fb | (3.43) 
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where f(x) is a Gaussian distribution function expressed as 
] x? ' 
i. > P| —— (3.44) 
./ 2x6 20 
and where 
N 
c= =e (3.45) 
Baer a — 4 : , a 
Using the transformation y = a and rearranging the terms in Eq. 3.43 vields 
0 N=1 
PLC) = | i: -2| 10 | gv)ay (3.46) 
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Using another variable transformation the inner integral becomes 
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Substituting this result into Eq. 3.46 gives 
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Figure 12. Diagram of Correlation Receiver for 2FSK/QPSK 
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= 
P,{C} = [i= 20(0+35 ) Poa (3.48) 
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Thus we have the probability of a correct symbol decision associated with the optimum 


receiver expressed in terms of the Q function. 


C. OPTIMUM RECEIVER BER 
In order to compute the Bit Error Rate (BER) of the optimum receiver shown in 
Figure 12, it is necessary to determine under what conditions a chosen bit will be cor- 
rectly indentified by the receiver for a given transmitted signal. The three bits assigned 
to each 2FSK/QPSK signal are identified as 1B, 2B and 3B corresponding to the most 
significant bit, the second-most significant bit, and the least significant bit respectively. 
We can assume without any loss in generality that s,(z) is transmitted. The probability 
that the 1B is correctly identified at the receiver, conditioned on the assumption that 

S(t) Was transmitted can be written as 
Pr{1B Correct | so} = Pr{r, = max rj|s} + Pr{ —r, = max 4159} 


Z 4 (3.49) 
+ Pr{r; = max 4,| 59} + Pr{ —rg = max r;| 59} 
j j 


Now 
Pr{r, = max rj| 59} 
j (3.50) 
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Mewesturther Condition on r,, that is r, = R, , since r,>—r, only for R,>0 , then this 


conditioning will be 


Pr{r,; = max rjl 59, r, =R, > 0} 
j 7 (3.51) 
= Pr{R, > kh, EG ye. R, > Fs, EO ee G06 ne ty Sos r, = R, > 0} 


Observe that 
Ri >Y%;, Rear => -R,<1r<R, (3i52) 


so that 
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Pee = max rls, t) = R, > 0} 


= Pr{-R,<r<R,, —Ry<1, Ry 2 oe oe Se _ 
Assuming s,(t) was transmitted, then from Eqs. 3.8 and 3.38 we have 
r(t) = JE $,(t) + n(d) (3.54) 
so that 
r=JE +n,, m)=m, =m, %=Ny (3.55) 


Due to the independence of the random variables n, , the set {r,} consists of independent 


Gaussian random variables having probability density functions (p.d.f.’s) 
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So that Eq. 3.53 simplifies to 
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Furthermore, from Eq. 3.57 
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Using this result with Eqs. 3.53 and 3.58 vields 
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so that integrating over all R, > 0 yields 
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Note that this expression approaches unity as = — oo as expected. Consider now 
Pr{ —r, = max r,| 59} 
j (3.62) 


= Pr{—r, >), —tp>—M, —tp >, —mh>r3, —Ig> 33, —tg >a ly, > —Tal 5p} 


If we further condition on r,, that is r,=R, , since —r,>r, can only be satisfied for 


R, <0, then this conditioning will be 


Pr{ —r, = max rjl 50, r, = R, < 0} 
/ 7 (3.63) 
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Let R’, = — R,, so that in terms of R’, , 
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so that 
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J 
2 ——— —_— 
R' R',+JE an Bs (3.68) 


Since r, takes on values R,, we can define a new random variable r’, = —r, where r’, takes 


@oeaies fk 4. 1.c,, i f;— A,, then r,— — R,= RF’, .. So the condition 
Yr, — aay Ke = 0 => =i 5 = ih 5, R’, > 0 = r’, = R'), R', > 0 (769) 


and given the symmetric p.d.f. for r, conditioned on 5, , we see that the conditional p.d_f. 


ofr’, is no different than the conditional p.d.f. of r, , so 


Pr{ —r) = max rj] 59} = Pr{r’, = max r;1 50} 
j / 


| ae ea poe OTE 3.70 
-| [1 -20( $2) ] fof 2%) -of BSE) fica | 
0 


Through a transformation of variables this expression is reduced to 


Pr{ —r, = max | 5} 
j 


i 


B 2 “ies oie 203TH 
=| P= 2008 [1-o(r+ a )-2-35)| = exp - Ja 


. ‘ : E ? 
Note that this last expression approaches zero as Pail Consider now 











Pr{r3 = max rj1 59} 
i) 
= Prity > 0), Fy > tp, 03> Tp, Ty > hg) Ty > 0g, Ty > tay Fy > 1g 15g} G72) 


= Pr{r, <rs, E> 83; rr, < 13, yy Be ha Or 0) rg <1, feS=ryls,} 


If we condition on r, that is r,=R, , we must have R,>0, or we cannot satisfy the 


2r; >0 condition. Thus 


Pr{r; = max 4,| 59, ry = R, > 0} 
J . (3.73) 
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From the mathematical form of this equation, it is obvious that this probability, when 
averaged over all R; > 0 produces the same probability as that given by Eq. 3.71, which 


also is identical to Pr{ —r, = max r,|s,} . Thus, we obtain from Eqs. 3.49, 3.61 and 3.71 
d 
Pr{1B Correct | so} 
fone) =< 3 
ve 
= l= 20 5 eae 


ve (3.74) 


o 


+ 3]"01 - 2000 | -o(v+ se )-ov- se ) from 


If the 2B is selected then 














Pr{2B Correct |S} = Pr{r, = max rjl so} + Pr{r, = max rl 59} 3.15) 
+ Pr{r; = max rj] 5} + Pr{r, = max 4jl 59} , 
j j 


The first term on the right hand side of Eq. 3.75 is given by Eq. 3.61, while the last three 
terms are all of the same form as Pr{r; = max r,|5,} from Eq. 3.72 and will all yield the 
same probability as given by Eq. 3.7]. TTanere Fone. the probability that the 2B 1s correct 
is the identical to the probabilitv that the 1B is correct. It can be shown in similar 
fashion that the probability of correctly identifving the 3B is also equal to that for the 
1B. Thus Eq. 3.74 is also the expression for P,{Correcr|s,} the average probability of 
decoding a bit correctly given that s,(z) was transmitted. 

Now consider the case of a given signal s,(t) being transmitted, where s,() can be 
any one of the signals s(7), i=0, 1, ..., 7. For any chosen bit. the probability that 
this bit is correct conditioned on the hypothesis that s,(1) was transmitted can be ex- 


pressed in a form 


fal 
Pr{ Chosen Bit Correct|3,} = Pr{r;= max r/13,} +), Pr{tm= max 115} (3.76) 
J nel ei 
msl 


where r, and the r,, terms are the signed correlator outputs +r,. The output r, corre- 
sponds to s, in that if r, is the maximum of all the +r, , then the receiver decides s,() 
was transmitted. The outputs r, correspond to three other s, which if selected by the 


receiver, would result in the chosen bit still being correct. The first term on the right 
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hand side will again vield the same probability as that given by Eq. 3.61. The terms in 
the summation are all of the same mathematical form as Eq. 3.72, and so are equal to 
the probability given by Eq. 3.71. Thus Eq. 3.74 is the expression for the probability 
that any chosen bit is correct for any given signal transmitted, and therefore vields the 
overall average probability of a correct bit, P,{C} . 

Some insight into the total probability of a correct bit decision P,{C}, may be gained 
by comparing Eq. 3.74 with the probability of a correct symbol decision P,{C} given by 
Eq. 3.48. Since N, the number of basis functions is equal to 4, the first term in the ex- 
pression for P,{C} is equal to P,{C}. This should come as no surprise, since the deriva- 
tion for the probability of the 1B being correct began with Eq. 3.43 where the first term 
is in fact the probababilty that this bit 1s correct because the receiver made the correct 
symbol decision. 

The second term in Eq. 3.74 represents the probability that the bit is correctly 
identified even though the receiver has made a symbol error. Recall from Figure 10 that 
only those signals which have vector representations orthogonal to each other have bits 
in common, and for 2FSK;QPSK there are three such signals meeting this requirement 
for any one fixed signal. The probabilities of erroneously choosing one of these signals 
over say S,(z) , are all equal as discussed above. Furthermore, the conditional probability 
of (incorrectly) choosing anv of the signals orthogonal to s,(r), can be shown to have the 
mathematical form of Eq. 3.72. Thus Eq. 3.71 specifies the probability of erroneously 
choosing anv one of the signals orthogonal to s,(z) . 

iitevabove results for 2FSK’QPSK can be extended to MFSK/QPSK, that is, a 
multi-frequency’four-phase signaling scheme. The structure of the optimum 
MFSK, QPSK receiver is the same as that of the general correlation receiver shown in 
Figure 11 on page 30 where the bank of N correlators corresponds to a pair of 
correlators, one each for the In-phase and Quadrature channels, for each of the K, fre- 
quencies. The expression for the conditional probability of a correct signal decision has 
already been given in Eq. 3.48. The expression for the conditional probability of decid- 
ing in favor of a signal that is orthogonal to the transmitted signal, vet does not result 
in a bit error takes on a form similar to that derived for the 2FSK/QPSK case. To 
demonstrate this, assume that s,(t) has been transmitted and examine the probability 
that one of the signals orthogonal to s,(r) is erroneously choosen over 5s,(r) . This corre- 


sponds to the probability that for some r,, A=2, 3, ..., N 
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Pr{r, = max r;| 59} 
j 


(3077) 
= Pr{ry, > Py Ty > Hl yy Np Toy tip py ey > See ee ee 
Further conditioning on r, yields 
Pri{r, = a jl 59, ip Ris R,> 0} 

aN. 

=| [Prt -Ry<1y < Rel %, r, = R;, R, > 0} (3.78) 
j=l 
J#k 


The first term in the summation (j= 1) is the same as that given by Eq. 3.66 for 
2FSK/QPSK. All the other terms (= 2, 3, ..., N; j#4&) are given by Eq. 3707s 
that 


Pr{r, = max rl 59, ea (OS R, > 0} 
j 
Ne? = = 
Real Rivest R',-VJE (3.79) 
= E Z 20(=)| c = “of Fa 
Now, integrating over all r, > 0 yields 


2 max r;| 59} 
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oo N-2 JE Wes 1? (ae 
[rsa fee B)-e ) tel BP 


sf oh 











Since there are \f — 2 signals orthogonal to s,(z) , and for any chosen bit, half of these 
signals have this chosen bit in common with the one in s,(), yields i = Naan 
ways of making a symbol error while correctly indentifying the bit in question. Thus the 
expression for the probability of a correct bit decision conditioned on the signal s,(2) 


being transmitted, is given by 


40 


P,{ Correct | so} 


i a 
= E =20ye FZ 1 gly)ay 
ae (3.81) 
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20 JE cE 
+0v~ Df C1200) 1- ofy+ z )-of aaa ) leon 








Using steps similar to those used in the analysis of the 2FSK/QPSK receiver as well as 
those presented above, it can be shown that Eq. 3.81 is also the expression for the overall 
average probability of a bit being correctly identified by the MFSK/QPSK optimum re- 
ceiver P,{C}. 

Since the integrals involved in calculating the bit error probability for the 
MFSK/QPSK optimum receiver cannot easily be expressed in closed form, the BER can 
onlv be evaluated using numerical techniques. Appendix B presents a computer pro- 
gram in FORTRAN language used to compute the actual bit error probabilitv for the 
optimum receiver by emploving numerical integration. An upper bound for the receiv- 


er’s bit error probability is given by Ref: 6 [pp. 255-257] as a union bound, namely 


ail kEs 2kE, 
Piles | (AL 2)Q No +O No 


This bound is a good approximation to the actual BER for & > 3 at moderate to large 
E, 

Ny 

D. DBD RECEIVER FOR MFSK/QPSK 

Another receiver structure for MFSK/QPSK signals is implementable by using DBD 








(3.82) 


techniques. Ref. 2 [pp. 19-25] A diagram of the structure of this receiver is shown in 
Figure 13 on page 42 for the 2FSK'QPSK case. The DBD receiver uses the same 
correlator subsystem as the optimum receiver Figure 12 on page 32, however the two 
receivers differ in the way the correlator outputs are processed. Whereas the optimum 
receiver processed these correlator outputs through a weighting network and then se- 
lected its largest output, the DBD receiver routes the correlator outputs through a two 
stage summing network and then compares the summer outputs to a Zero threshold in 


order to arrive at binary decisions. 
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E. MFSK/QPSK DBD RECEIVER BER 

The analysis for the derivation of the BER performance for the 2FSK/QPSK DBD 
receiver is similar to that carried out for the performance analysis of the optimum 
(symbol discrimination) receiver. We begin by assuming s,(2) is transmitted. Since the 
correlator outputs are identical for both receivers in question, the analysis on the per- 
formance of the DBD receiver can proceed with the correlator outputs, which are again 
characterized as independent Gaussian random variables with density functions given 
by Eqs. 3.56 and 3.57. Note from Figure 13 that there is a path from each correlator 
output to the input of all second stage summers. The outputs of these summers, denoted 


peepee? = 1, 2, 3 , are given by 


b, =r, -—r, +h by 
bo =r, +r,+134+1%%, (3283) 
Dy =r — 9 + yt Fy 


Each b, will also be a Gaussian random variable having a mean equal to the sum of the 
means of the r,, and due to the independence of the r,, the variance of b, will consist of 


N 
the sum of the variances of the r, , each of which equals =e Thus we have 


b; ~ AGE, 2N,) 
b, ~ N(/E, 2No) 3.84) 
b,~ Al —/E . 2N,) 


where Aa, 5) is used to denote a Guassian probability density of mean a@ and variance 
6. With the statistics of the b, specified comparing these random variables to a zero 
threshold yields the following conditional probabilities of decoding the bits correctly 


namely 


(x= a. 


Pr{1B Correct|s,} = Pr{b, > 0} -| == | — aN) las 
0 ° (3.85) 


af 2n(2Np) 


d-JE)--(JE) 
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Pride Correct| so} = Pr{b, > 0} = ie exp| — Ges dx 
x 2x(2No) 2(2No) 3.86) 


E, 
vol of 2No ) 





(x+ JE) 
Pr{3B Correct | so} Be Pr{b, < 0} = . nar | = OI Je me 


Since each of the individual bits have identical probability of correct decoding, the cor- 





responding average probability conditioned on the assumption that s,(7) was transmitted 


Pa( Correct} = 1 ~ of = (3.88) 


and consequently the conditional bit error probability is 


1S 


P,{ Error |S} = 1 — P,{Correct| so} = of a (3.89) 
No 
Consider next the-case in which any one of the s(7}; i=O;4, ..., 71s tramsmittediion 


any of the receiver correlators. If the transmitted signal is matched to the local signal 


correlator, its output r, will be 
r= JE +n, (3.90) 
and for the correlator having antipodal local signal matching, the output will be 


r=—-JE +n, (3.91) 


All other correlators will have an output given by 


r=n; (3-92) 


Examining now any one of the second stage summer outputs b, , if the bit corresponding 


to b, for the signal sent is a logical zero, then b, is statistically characterized by 
baw NG/E ; 2N5) (3.93) 


In this case the bit will be correct if b, is greater than the zero threshold, and the prob- 


ability of the bit corresponding to b, being correctly decoded is given by 


ry vd, = VS === XK = Son meee 
| ./2n(2Np) 2(2No) 





(3.94) 
=o) ) 
However, in the case where the bit is a logical one, then 
b, we jE, 2Ne) (3.95) 


Now the bit will be correct if b, is less than the zero threshold, and the probability of 


correct decoding is given by 


0 ior 5) 
1 (x+ VE) 
Pr{b, <0} = ON “la 
. I s/ 2n(2No) | 2(2\o) | 


ee E 
a o( a 


For both cases the probability of a bit error is the same, and it is identical for all b, re- 


(3.96) 





gardless of which signal sz) is transmitted. Thus, we can write the total bit error prob- 
abilitv for the 2FSK/QPSK DBD receiver as 


P,{e} = of. oe ) ; (3.97) 


Extending the 2FSK/QPSK results to the general case of MFSK/QPSK in which 
there are M signals, N receiver correlators, k bits, and K, frequencies where the 


biorthogonal signal relationships of Eqs. 3.35 - 3.37 again hold, for each of the random 
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variables b, , n=1, 2,..., k , the corresponding variance will be the sum of the variances 


of the n,, so 


N 
V a No an NNo = . 
ar{b,} = es om KyNo (3198) 


j= 
and again the mean of the b, will be +./E . Following the same steps as those carried 
out in obtaining Eqs. 3.90 - 3.97, the total bit error probability for MFSK/QPSK DBD 


receiver 1s given by 
[2k 
Ee RTS 3.99 


F. BANDWIDTH OF MFSK/QPSK 

The bandwidth of MFSK/QPSK signals can be derived from knowledge of QPSK 
and MFSK bandwidth requirements. Each carrier is switched at the svmbol rate R, 
\Wihere re . This gives each carrier a two-sided bandwidth of 2R, Hz. To maintain 
orthogonality between signals at different frequencies in MFSK modulation, the fre- 
quencies are separated by R, Hz. This same relationship must also be maintained in 
MFSK/QPSK. The spectrum of the 2FSK;QPSK signal is plotted in Figure 14 on page 
47. From the graph, the total null-to-null bandwidth of the 2FSK/QPSK signal can be 


seen to be 
W = 3R, (3.100) 


Proceeding to the family of MFSK/QPSK signals, Figure 15 on page 48 shows the 
spectrum for such signals. By examining this figure, the total null-to-null bandwidth can 
be shown to be 

W = (Ky+ ))R; 


Seal (3.101) 


This general result for MFSK/QPSK correctly yields the result for the special case of 
2FSK/QPSK in which k = 3, namely 


W = (2°? + 1)R, =3R, (3.102) 


46 


*‘syusuodulos {eloads 


JeNprAIpUl MmoYs saut] payseq :YSdO/MSAT 10J winaysadg apnjziuseyy 


AONANGAMA 





‘pl aansly 


00 


20 


yO 


9°0 


g°0 


Oy 


HONLINOVAW GaZrIVNUON 


47 


‘sjuauoduios jeijoads 


[ENPIAtpur MoYs sour] payseq :WSdO/MSAIN 40J umapadg apnyuseyy 
AONANGOAYA 





ST eansty 


20 0°0 


v0 


9°0 


#°0 


Ad NLINOVA GHaZTTVWUON 


aT 


48 


which agrees with Eq. 3.100. 


The bandwidth efficiency of a digital modulation scheme is given by Ref. 6 [pp. 
393-395] as 


Ro 
aT (331103) 


and is one way to evaluate how efficiently different modulation schemes use the available 
channel bandwidth. 


For M-ary PSK, bandwidth efficiency is given by Ref. 3 [pp. 
211-213] 





3.104 
R (3.104) 


and for M-ary FSK the bandwidth efficiency is 


Ry 
TESK™ (M+ DR, 
kR, : 
UitpR NR, (GAl05) 
k 
i hey 





For MFSK/QPSK, using Eq. 3.101 we obtain a bandwidth efficiency given by 


Rp 
("al WR. 
KR, 

eae 
ee: 
. aa 


MFSKIOPSK = 


(3.106) 


A table of bandwidth efficiencies for M-ary PSK, M-arv FSK, and MFSK/QPSK is 
given in Chapter 5. 
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IV. RAYLEIGH FADING 


A. THE RAYLEIGH FADING MODEL 

Often, when a signal is transmitted through a channel in which multipath propa- 
gation and/or heavy precipitation conditions exist, fading occurs with a consequent 
change in the signal amplitude. One way to model this fading is to redefine the signal 
amplitude 4 which up until now has been treated as a constant, as a Rayleigh distributed 
random variable A with p.d-f. 





A A? 
ii(A) = exe u(A) (4.1) 
5 Ae aye 
where u(-) is the unit step function, and from Ref. 7 [pp. 110-112], the value of A, is re- 


lated to the second moment of A as 
E{A’} =2A2 (4.2) 


In order to evaluate the effect of Ravleigh fading on the BER for a given modulation 
scheme and corresponding receiver, the expression for BER under nonfading channel 
conditions must be averaged over the distribution of the signal amplitude. This requires 
expressing the bit error probabilitv as a function of signal amplitude, weighting the 


function with the amplitude p.d.f. and integrating this product over all values of A. 


B. EFFECT ON BER OF M-ARY PSK USING DBD RECEIVER 
For the 8-PSK DBD receiver, the (non-fading) bit error probabilities are given by 
Eqs. 2.1 - 2.2. These probabilities must now be expressed as functions of the random 


aniplitude A. From Eq. 1.1 it is apparent that the signal amplitude is given by 


A= Je (4.3) 


so that 





E= , (4.4) 


Since the signal amplitude is now a random variable, the signal energy also becomes a 


random variable E, which is expressed as 
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p= (4.5) 


Substituting this result into Eqs. 2.1 through 2.3 vields 


Pr{1B incorrect} = 
ae | Ts ill Ts (4.6) 
5 (1 N, COs >) ey of Ny sin 2) 
Pr{2B incorrect} = 
Te Te ; 
QUA Ty [cosa+sina]]+Q\A TN [ cos a — sin a] 
2 YO (4.7) 
i ° T, . 
—2Q0\ A —[cosa+sing}]Q( A —- [cosa —sin a] 
2N5 2No 


and as before, 














Pr{3B incorrect} = Pr{1B incorrect} (4.8) 


The bit error probability under fading conditions for the most significant bit, namely 1B, 
is obtained by weighting the probability given in Eq. 4.6 by the p.d.f. of the random 


amplitude A and integrating over all A. That is 


Pray{ 1B incorrect} -| Pr{1B incorrect}f,(A)dA 


a _ ae (4.9) 
2 Fos iF cosa) 440( Ne ins f,\(A)dA 


This integral can be solved using the identity given by Ref. 8 [pp. 12-13] as 








| POC (Aad = 4) 1-1 (4.10) 





Using Eq. 4.2 and the linearity of the expectation operator we can write an expression 


for the average signal energy under fading conditions, namely E{E} = E as 


5) | 


2 2 
2 


E{A?} = 3 (243) 


= A 
E=£ 





(4.11) 


v5 


te 


Using Eq. 4.10 in order to evaluate Eq. 4.9 vields the bit error probability under Rayleigh 


fading, namely 
Pragy{1B incorrect} = 


| 2 eee eee 


eee 

2 4 = at = 
(  cos*s) al 
No 


Next, the BER under fading conditions for the second-most significant bit, namely 2B, 


(4.12) 





is obtained by integrating the product of the BER given by Eq. 4.7 with the Rayleigh 
p.d.f.. The contribution from the first two terms of Eq. 4.7 can be derived using the 


same steps as those carried out above for the evaluation of P,,,{1B incorrect}. Since the 


Rey 
third term in Eq. 4.7 involves the product of two Q functions, a new identity is required. 


That is 


| ocnoayy(aaa = 
| ] co 
ie aa ” a + a é 
Bran”! + a tan 7 a (4.13) 
] 


+7 a GGG 
2n P+ 2n e+ 
V A V Ae 


the proof of which is presented in Appendix A. Using this identity for the contribution 
of the third term of Eq. 4.7 and combining this with the contributions from the first two 


terms of Eq. 4.7 yields 


Nn 
tw 


2 {2B pe = 


azar 


ear [mae (4.14) 
u, tan” || ————_ 








a) 


om 


-1 
cane Ui; 


 afae{Ey  afafEy 


where u, and uw, are given by 


u,; =cosat sing 


(4.15) 
Uo = COS a — SIN & 
The overall average bit error rate under fading conditions can be obtained from 
Py, Raye} = = Proy{1B incorrect} + = > Pray(2B incorrect} (4.16) 


For the 16-PSK DBD receiver operating under Rayleigh fading conditions, much 
of the performance analysis can be carried out using the same procedure as the one used 
for the evaluation of the performance of the 8-PSK DBD receiver operating under simi- 
lar Rayleigh fading conditions. The bit error probabilities under non-fading conditions 
are given by Eqs. 2.4 - 2.6. Rewriting these equations in terms of the random variable 


A, for the most significant bit we obtain 





Pr{ iB. imeorrect} = (4.17) 


|— 
to 
a 
> 
2|n 
& 
eae 


where the variables v, are defined as 


¥, = COS & 


r, = cos £ i 
vz = sing (4.18) 
v, = sin B 


and for the second-most significant bit 


Pr{2B incorect} = 


FE) a) (Va) ay 
-o(sl te) 


where the variables w, are given bv 








~ 








w, = cosa — sine 
Wy = COS & + Sin & 
w; = cos 6 — sin Bf 
Ww, = cos B + sin B 


(4.20) 


For the third-most significant bit 


Pr{3B incorrect} = 


i= / ze 2 
i A = i ool - = ils a, ear a, B) F 
i=) hie er a3 
Te 
= of. [x, alt a phar “em 
2 oe 2 
Te ca . A & 2, 
+ DAY 2nNo I exo( - (OBAa ear a, B)e 
V No 


where the functions f(T, «, 6) and g(I, a, £) are defined in Table 1 of Chapter 2. For 
the least significant bit, from Eq. 2.3 it is obvious that 


| 
Ts 








rf 


ro |— 





rf— 


Pr{4B incorrect} = Pr{1B incorrect} (4.22) 
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For the most significant bit, the BER under fading conditions is obtained using Eq. 4.17 


and the identity given by Eq. 4.9, to yield 


(4.23) 





where the variables v, are defined in Eq. 4.18, and the last equality in Eq. 4.23 follows 
from Eq. 4.22. The BER under Ravleigh fading conditions for the second-most signif- 


icant bit is obtained using the identities of Eqs. 4.9 and 4.13 to vield 


ew 





a! 
Wa w, tan 


Lr Sil rE = 
2nd +2{ an, vd +2{ 
i\o 2G) 


where the w, terms are defined in Eq. 4.20. For the third-most significant bit, the inte- 


grals to be evaluated are all of the basic form 


So) 





| 4] : eo( - 
—oo *0 
"7 1 % 2 2 2 
=| gf, a, B) ya Ih ee exp( — x°)dxaT (4.25) 
c 
0 


E; Ag 


7 1 Ja 
-| eis a, B) avy 4 |« 


Jeu a, pal f,(A)adA 











0 


where ¢, is defined as 





(4.26) 


and 





if 
P= Ty IAT, a 8) (4.27) 


The variable x is defined as x = Ae,, so that the evaluation of the integral of the product 


of the first term in the summation of Eq. 4.21 with the Rayleigh p.d.f. becomes 


[2 si [Pex [Fexo( - 5s 2 ae a, B)dT fy(A)ad 
1 (tas ae ee 
ay ic = oS [ian =e av (4.28) 


je = gl, ~ 9 
il E| lst BY ip 
cl Ee FAT, «, neil? 
No 


Multiplying the second term within the summation of Eq. 4.21 by the Rayleigh p.d.f. and 
integrating, yields 








q 














[4 ze = ‘exn( -4 oi Zs AT, a, pet. O, ola / Cen ~ BT a, faa 
on ie Ae 
=|" A 22No i ool - 2No 
lg g(x)dxdTf,(A)dA (4.29) 
4% aT, a, B) 


s A , A ye 
= é,exp( — A’p; a exp Piaaxd 
Pat ia ¢; exp( pig (x)u(x — Ag;) re = ve )s )dAdx 


Where €, is defined as 
Ie 
¢; = x ol, a, B) (4.30) 
‘No 


Carrying out the integration with respect to 4 first, while realizing that u(x — A&)u(A) 











Jean a, B) a 








é : AG : 
is nonzero only for A in the rangeO<A< 2 vields 


al 


Xe 5] 

Peal A 

ea exp( — A 29) A exp dA 
\ \/2n “Vie eo 


0 0 
= 2 
= |’ == exp( — aA dA - exp| — sii 
Ag! 2z F 2a; Dak e 


where a, is defined as 














a= pt > 


ase 
ve (4.32) 


so that Eq. 4.29 becomes 


S7 


1 [*[" ew -4). 
Asean 0 “0 ' /2n - 
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x 2 
Si a; 
| els exp( — a,A°)dA -- TaF, = = =) dxal 
0 i 








24; isi 
| : (4.33) 
é, 252) 
. —— eo(-)| so exp( — aA’)dA 
Deal sf] V22 0 dxaT 
ee 2 : : 
ON aco | |) Se ox 
2a f/2n 2 gE 


Focusing on the integration with respect to x only, we have for the first integral with 


RESPEC LOX, 





x 
co 2 Nee 
| — exp( = =| re exp( — a,A°)\dAdx 
a wanes i 


a os x? 1 2 ala 2 
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and for the second integral with respect to x 





—— exp( = Ls )és = = (4.35) 
|, af 2n 2 2ay;/2n 


where y, is defined as 
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Thus, Eq 4.29 becomes 


“Ay 7 fe *axo( - Fe aN, 
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From our previous definitions we can rewrite the expressions for a, and y, as 
a;=/ + 
i t te 
(4.38) 
as 
me 7 5 
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Observe that the terms in Eq. 4.37 can be expressed in terms of the ratio of average 


signal energy to noise power spectral density as 








z E g; | ie ’ 
g; ae 1. = e gi, af) : (4.40) 
4Agq2 | eva. 0, pif? 
No 
a 
a (4.41) 
J 2a; 





59 


and 


- mills oe (Ae 
a (4.42) 


Se 5 Qa, B) 
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Using these results, the first term on the right hand side of Eq. 4.37 becomes 


l * ci eee 
Qn i 2 ae No an z (4.43) 
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While the second term on the right hand side of Eq. 4.37 takes on the form 
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Finally, from the third term on the right hand side of Eq. 4.37 takes on the form 
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so that the bit error probability for the third-most significant bit 1s obtained from the 
sum of Eqs. 4.28 and 4.43 - 4.45 . These integrals have not been expressed in closed 
form, however corresponding numerical evaluations can and have been carried out 


without difficulty. The total BER under Rayleigh fading conditions is finally given by 


Py, Raye} = > Pray{1B incorrect} + oa Pray{2B incorrect} 
! (4.46) 
7 Pray{3B incorrect} 


Appendix C presents a listing of the FORTRAN computer program used to evaluate the 
BER under Rayleigh fading conditions for the 16-PSK DBD receiver. Data from this 


program is plotted and presented in Chapter 5. 


C. EFFECT ON BER OF MFSK/QPSK USING DBD RECEIVER 
Recall that for the MFSK;QPSK DBD receiver, the BER is given by Eq. 3.99 as 


P,{e} = of — (4.47) 


This expression must be rewritten as a function of the amplitude A in order to determine 
the effect of Ravleigh fading on receiver performance. Substituting Eq. 4.4 into Eq. 4.47 


vields 


(4.48) 





Weighting this probability by the Rayleigh p.d.f. and integrating over the variable A, the 


average bit error probability in the presence of signal fading is 


oo de 
a Ray(€} -| of 4, NNo tana (4.49) 


This integral can be evaluated using the identity of Eq. 4.9 as 
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l _ 
Py, Ray(?} — is (4.50) 








D. EFFECT ON BER OF MFSK/QPSK USING THE OPTIMUM RECEIVER 
Evaluation of the (non-fading) BER for MFSK/QPSK using the optimum receiver 
has been carried out in Chapter 3, Section C. Yet a different approach to determining 
the BER under Rayeigh Gane conditions becomes more practical than using modifica- 
tions to the approach previously used to evaluate the performance of the DBD receiver. 
Since the BER given in Eq. 3.81 as a function of s, the ratio of signal energy to noise 


variance, namely 














s=-—5 (4.51) 
where 
at 4,52 
Os a ( : ) 
has been previously evaluated, a transformation of variables is suggested to yield 
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where the range of sisO<$<oo. Since the transformation is linear, the p.d.f. of s can 


be written as 














oT (4.54) 











where S? is defined by 





The BER under fading conditions can be evaluated by averaging the BER given by Eq. 


3.81 over the range of the new random variable s as 
Te Raye) oa 
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This integration also requires numerical techniques, however the program presented in 

Appendix B can be used as a subroutine for the program listed in Appendix D which is 

used to evaluate Eq. 4.56. Data from the program of Appendix D is plotted and pre- 

sented in Chapter 5. 
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V. RESULTS 


A. PERFORMANCE COMPARISONS 

The symbol error rate (SER) of the receiver for M-ary PSK signals is given by Eq. 
1.4, while Eq. 1.8 is used to evaluate the SER of the receiver for M-ary FSK signals. 
The SER for MFSK/QPSK signals when using the optimum receiver of Figure 12 on 
page 32 is given by Eq. 3. 48. Table 2 lists the ratio of bit energy to noise power spectral 
density (PSD) level, = : 
ary PSK, M-ary FSK, and MFSK/QPSK signaling schemes with k equal to 3 and 4 
bits/symbol]. The table shows the MFSK'QPSK signaling scheme yields slightly better 
performance than M-ary FSK with a nominal 0.2 dB difference, in terms of = és , while 
outperforming M-ary PSK by an average = difference of 52 cRior.—s a 10.4 dB 
fork =4. 


required to achieve symbol error rates of 10-3 and 10-5 for M- 


Table 2. PERFORMANCE COMPARISONS AT KEY SYMBOL ERROR 
RATES 







EB, 
Bits per =e (dB) Required to Achieve 
Symbol “0 
oe 10> SER —-10-$. SER 


[_iersk Ss“) «dD 6 














Modulation 





The BER for MFSK/QPSK signaling using the optimum receiver shown in 
Figure 1] on page 30 is evaluated using Eq. 3.81. Appendix B presents the FORTRAN 
computer program used to compute this optimum receiver's bit error probability by 
employing numerical integration in order to evaluate Eq. 3.81. That is, P,{e} along with 
the union bound approximation to this probability as given by Eq. 3.82 is computed for 


varying values of ak The program is also used to evaluate P,{e} for the MFSK,QPSK 
ia | 


DBD receiver, as given by Eq. 3.99. Data generated by this program is plotted in 
Figure 16 for k =3, Figure 17 fork =4, Figure 18 for A=5, and Figure 19 fork =6. 

In all of the above referenced figures, numerical values for the BER obtained using 
the union bound approximation yield plots eae coinciding with the actual BER plots 


for the optimum receiver at moderate to high —— values. This fact helps to confirm the 


appropriateness of the results of the ae integration used to evaluate the actual 
receiver BER. From these plots, it becomes evident that the MFSK/QPSK DBD re- 
ceiver suffers a Pe ae penalty when compared to the performance of the optimum 
receiver. Table 3 lists x ratios required to achieve BER’s of 10-3 and 10~ for both the 
optimum and the DBD receiver structures for MFSK/QPSK modulation. The optimum 
receiver has an average 2.5 dB better performance than the DBD receiver for_k =3 
bits‘symbol, which increases to an average value of 5 dB for k =4 in terms of dif- 
ference. For k = 6, the optimum receiver has an average 10 dB advantage over the DBD 


receiver. 


Table 3. PERFORMANCE COMPARISONS OF MFSK/QPSK RECEIVERS 





Modulation’ Bits per 
Receivers Symbol 
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The reasons for this increased performance penalty can be explained as follows. 
First, as can be seen bv examining the DBD receiver BER curves as a function of oe 


. . z 0 
the performance of the DBD receiver degrades as k, the number of bits per symbol is 
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Figure 16. 2FSK/QPSK Bit Error Performance (k = 3 bits/symbol) 


66 


10° 


1 
S 
— 


BIT ERROR PROBABILITY 


on 


LEGEND 
_ OPTIMUM RCVR 
DBD RCVR 





~10 — 0 10 - 15 


5 
EB/NO (DB) 
Figure 17. 4FSK/QPSK Bit Error Performance (k = 4 bits/symbol) 
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Figure 18. S8FSK/QPSK Bit Error Performance (k = 5 bits/symbol) 
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increased. Recall from the BER analysis for the DBD receiver, that the variance for the 
summation output b, is the sum of the variances of the correlator outputs r,, When the 
number of bits ( and consequenctly the number of correlators ) is increased, there is a 
linear increase in the variance of b,. This variance increase, which is strictly noise de- 
pendent, causes the bit error probability increase also. On the other hand, examination 
of the optimum receiver bit error curves reveals that, with increasing k, the optimum 
receiver exhibits improved performance. This improvement reflects the tradeoff that 
exists between bit error -rate performance and signal bandwidth, as increasing k results 
in a greater signal bandwidth. Note however that this improved receiver performance 
would not be apparent if the bit error probability were plotted against unnormalized 
SNR instead of bit-normalized SNR, that is, x, Ref. 6 [pp. 172-175]. 


B. BANDWIDTH EFFICIENCY 

The bandwidth efficiency for M-arv PSK, M-ary FSK, and MFSK/QPSK signaling 
is given by Eqs. 3.105, 3.106 and 3.107, respectively. Bandwidth efficiencies given in 
bits’ Hz for values of k ranging from 2 to 6 bits/symbol are listed in Table 4 for these 
three signaling schemes. The table shows that M-ary PSK is more bandwidth efficient 
than both the M-FSK and MFSK’QPSK schemes. MFSK/QPSK is from three to four 
umes as efficient as M-ary FSK, with the greater differences occuring at higher values 
of k. MFSK/QPSK operates at 67% the efficiency of M-ary PSK for k = 3 bits’symbol, 
and at about 40% the efficiency of M-ary PSK for A = 4 bits’symbol. 


Table 4. BANDWIDTH EFFICIENCIES FOR M-ARY PSK, M-ARY FSK, AND 
MFSK/QPSK 


Number Number Bandwidth Efficiency 
of a of Signals MPSK MFSK QPSK 





C. EFFECTS OF RAYLEIGH FADING 
The BER of the 8-PSK DBD receiver operating under Rayleigh fading signal con- 
ditions is given by Eq. 4.16 as a function of average =o . This equation has been used 
0 
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to evaluate receiver performance as a function of average = for various values of the 
angle alpha. The plots shown in Figure 20 on page 72 demonstrate that the best re- 
ceiver performance occurs for o = 22.5°. The curves for o = 20.0° and « = 25.0° and for 
a = 17.5° and o = 27.5° are coincident, and therefore appear as one line for each pair. 

Appendix C presents the FORTRAN program used to evaluate via numerical inte- 
gration the bit error probability under Rayleigh fading signal conditions for the 16-PSK 
DBD receiver. Data generated by this program is plotted in Figure 21 on page 73 for 
various values of the angle « while, maintaining the relationship B = 45° —a. This plot 
shows the best receiver performance is achieved for a = 11.25° . Curves for a = 9.75° 
and a = 12.75° as well as for a = 8.25° and a = 14.25° are coincident and appear as only 
one line. 

The receiver's BER for MFSK/QPSK using the DBD scheme, is given by Eq. 4.50 
when Rayleigh fading signal conditions exist. Figure 22 on page 74 illustrates the plot 
of this BER for values of k ranging from 3 to 5S bits/symbol. The best performance is 
obtained for k equal to 3, and performance degrades with increasing values of k. The 
performance curve has a near linear characteristic for values of average 7 exceeding 
about 5 dB as opposed to the “waterfall” characteristic of the BER curves when Ravleigh 
fading conditions were not present. 

The BER for MFSK/QPSK signaling using the optimum receiver under Ravleigh 
fading signal conditions is given by Eq. 4.56. Appendix D presents the listing of the 
FORTRAN program used to evaluate this BER, which involves double numerical inte- 
gration. Data generated by this program is plotted in Figure 23 on page 75 for values 
of k ranging from 3 to 5 bits’svmbol. The curves show a somewhat surprising result in 
that the best performance is obtained when k is equal to 3, and performance degrades 
with increasing values of k. This result is opposite to that encountered when Rayleigh 
fading conditions Were not present, in which increasing k produced improved receiver 
performance. The curves of Figure 23 also exhibit nearly linear characteristics for val- 
ues of average x greater than approximately 0 dB. E 

Table 5 on page 76 presents a summary of the average x values required to 
achieve bit error rates of 10-! and 10~ for all of the receivers discussed above, under both 
non-fading and Rayleigh fading signal conditions. For the MFSK/QPSK receivers, in 
order to simplifiv the comparisons, only the results for k=3 are presented. For the 
MFSK/QPSK receivers, the performance degradation is about the same for both the 
optimum and the DBD receivers. In order to achieve a ZENS of 10 both receivers sulfer 


a nominal 3 dB performance loss in terms of average —— = under Ravleigh fading condi- 
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Figure 20. 8-PSK DBD Receiver Bit Error Performance in Rayleigh Fading 
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Figure 21. 16-PSK DBD Receiver Bit Error Performance in Rayleigh Fading 
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Figure 22. MFSK/QPSK DBD Revr. Bit Error Performance in Rayleigh Fading 
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tions, and performance is correspondingly degraded by about 10 dB at a BER of 10°. 
For the M-ary PSK DBD receivers, performance is degraded by about 4 dB and by 12 
dB at a BER of 10-' and 10-?, respectively. 


Table 5. EFFECT OF RAYLEIGH FADING USING VARIOUS RECEIVERS 

































Modulation; Rayleigh ra (dB) Required to Achieve 
Receiver Fading 10” BER 10-? BER 
[__. _2Fskorsk |__|) = 
Optimum | No | 6 | 
Receiver 
[___rskgpsk |) 
Se a 







Receive 
S-PSK | ae 

5 _ Ne |__| aaa 

Receiver 

1G-PSK a 

oa 

Reveiver 










76 


VI. CONCLUSIONS 


The results presented show that MFSK/QPSK modulation is superior to M-ary FSK 
in aspects of error performance and bandwidth efficiency. While MFSK/QPSK has only 
slightly better error performance over M-ary FSK, MFSK/QPSK requires only one-third 
to one-fourth the bandwidth used by M-ary FSK signaling. Also, since the 
MFSK/QPSK receiver requires the use of only one fourth the number of synchronized 
carriers and integrate and dump circuits required by the M-arv FSK receiver, it can be 
argued that the MFSK/QPSK receiver is less complex than the M-ary FSK receiver. 
The results also show that MFSK/QPSK is comparable to M-ary PSK in the sense that 
it offers improved error performance over M-ary PSK, but does not reduce bandwidth 
efficiency as much as orthogonal signaling. 

In comparing the different types of MFSK/QPSK receivers, the DBD receiver offers 
the ability to decode bits individually without symbol-to-bit mapping circuitry; however, 
this is accomplished at the expense of error performance. This error performance pen- 
alty becomes prohibitive for larger values of k. Since the advantage of receiver simplicity 
is no longer as attractive for larger k. the DBD receiver should be considered as an al- 
ternative to the optimum receiver only for smaller values of k. 

Analysis of the effect of Rayleigh fading in M-ary FSK DBD receivers and in 
MFSK/QPSK receivers has generally demonstrated severe performance degradation in 
the sense of producing large increases in the bit energy to noise PSD ratios required to 
maintain only modest BER performance. The penalties increase rapidly for lower bit 
error probabilities. Unless fading effects can be overcome, all of the receivers studied 
become impractical in the presence of fading. The MFSK/QPSK receivers show a slight 
reduction in performance degradation over the M-ary PSK DBD receivers. For the 
M-ary PSK receivers, the same value of the angle adjustment « which produced the 
minimum BER for a non-fading channel also produced the best overall performance 
under Rayleigh fading conditions. 
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APPENDIX = A. INTEGRAL IDENTITY 


This appendix presents the proof for the integral identity used in the text as Eq. 4.13. 





| onowayn(aaa = 
a 
[p+ > mer 
Btan™! eo o tan” 2 - (By 
2 ee 
4 
22 B+ l. 27 +t 
Ag Ao 
We integrate by parts, 1.e., 
| udv = uv|>° + | vdu (B.2) 
0 0 
with the following substitutions 
u= OaA)O(BA) (B.3) 
dv = f,(A)dA 
2 
= a exo = 5 yea (B.4) 
Ao 2A6 
Differentiating Eq. B.3 with respect to A, yields 
d : 5 
FT = Q(eA)OABA) + O(aA)O(BA) (B.5) 


where 
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ss i exp a dy 
dA ; Jin ‘ 2 (B.6) 








where the last equality is obtained using Liebnitz rule to differentiate the integral. Sim- 
ilar results hold for O(f4A), thus yielding 


Dee: 
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Now, integrating Eq. B.4 to obtain r results in 
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v= A exp = A aA 
Ao 2A9 











; (B.8) 
2A5 
Solving for the first term on the right hand side of Eq. B.2 yields 
A? . 2 
urls? =| — O(d)O(BA) exp| A> }| = 0700) 
2A / |o (B.9) 
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For the second term on the right hand side of Eq. B.2 , we have 


fore) foro) 2 
| vdu =| Expl — a = 10 
0 0 2A 


2 pes A Ak 
fo = exp( - pe ) +16) = exp( aS Ie (B.10) 


Nee 























Vet 


-| oa = exp( — 274°) + Q(BA) = oo 20°) a 
0 vee 


79 


where z, and z, are defined as 








oa 
2, = 5 ar Dae 
and 
2 
= jee 
22> 5 + 2A2 


The form of these integrals is similar to the identity given by [Ref. 9: pp. 1-20] a 
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where a and b may be complex, and 
ee Real(a) 
Imag(a) 


Using this identity and the equality 
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O®) = + Erfe( = ) 


produces an equivalent form, namely 





co 22 1 act+ b/2 
| O@xlexp@iw)aa— a rs 
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Using Eq. B.16 to solve Eq. B.10 yields 


0 B a+ jz /2 a B + jz /2 
vdu = ——= In} ——*— | + —4—= in| ——*S 
0 janz/2 ey 2 janz./2 B —jzyV/2 


(B.11) 


(B.12) 


(B.13) 


(B.14) 


(B.16) 


(B.17) 


Focusing on the argument of the first log function and converting to complex exponen- 


tial notation produces 


80 


a a: z,/2 
om vj a? + a eo tan”*{ aa ) 
a + j21/2 


a — jz,./2 Ape 
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where the last equality holds by restricting both o and z, to be positive, which is indeed 
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the case. Using this result in Eq. B.17 yields 
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Finally, substituting Eqs. B.9 and B.19 into Eq. B.2 yields the right hand side of Eq. B.1, 


thus proving the identity. 
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APPEN 


DIX  B. MFSK/QPSK BER COMPUTER PROGRAM 


Tete HeFe Fevewe we Fe Ie Fe Ie He Ie Me Ie Te He Ie Fe Ie Fe eT ETE Te Fe Te IE Te ENE TEE Me Fe Te FEN Te TF He Fe Te Te Te HIE eH CT Fe He ET TI I TT IIT 


* PROGRAM TO COMPUTE THE BIT ERROR RATES FOR THE BI-ORTHOGONAL 
oo SIGNAL SET (BOSS) FOR THE OPTIMUM RCVR AND FOR THE MFSK/QPSK 
“ DIRECT BIT DETECTION (DBD) RCVR. 


we 


dededesesede steve ese re dete Fede He Fees Te eFC Te eT We Ie eT EK Tee I eT I ee ese Fe He I ete Te Tek Kee ICI KK KR IK KI 


PROGRAM PRBTER 


IMPLICIT CHARACTER*1 (A-Z) 


EBNODB, INIT, FINAL, DLTADB 
EBDVNO,SREDSI 
PRBTCR, PRBERR, INTGL1 , INTGL2 , INTGL3 
DBDPBE , UBAPBE 
DCADRE ,F1,A1,B1,AERR,RERR,ERROR 
F2 AZ 46263 
PIE ,SQR2PI ,SQR8PI ,ONE 
Te Ro Ny tls Patek NOMEaS 
LEVEL, LEVOLD 


Lying oe: 
SREDSI,N 


INITIALIZE CONSTANTS 


(FINAL - INIT)/(NUMPTS-1) 


SUPPRESS ALL ERROR MESSAGES. 


DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
INTEGER 
INTEGER 
EXTERNAL 
COMMON 
C 
PIE = 3. 141592654D0 
SQR2PI = DSQRT(2*PIE) 
SQR8PI = DSQRT(8*PIE) 
K = 6 
M =e Ke 
N = 2**(K-1) 
F = 22 CK=2) 
A2 = 00. ODO 
Bl = 20. 0D0 
B2 = 20. 0D0 
ONE = 01. 0D0 
NUMPTS = 101 
TINGE = -10. 0D0 
FINAL = 8. ODO 
DLTADB = 
AERR = 0.0 
RERR = .0000001 
C 
LEVEL = 0 
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CALL UERSET( LEVEL, LEVOLD ) 
CALL XUFLOW(0) 


C COMPUTE PROBABILITY OF BIT ERROR 
C FOR VALUES OF EB/NO IN THE RANGE OF: 
Cc INIT < EBNODB < FINAL 

WRITE(8,100) K,F 

WRITE(8,*)' ' 

WRITE(8,*)' OPTIMUM RCVR DBD RCVR 

& UNION BOUND’ 

WRITE(8,*)' BIT ERROR BIT ERROR 

& APPROXIMATION' 

WRITE(8,*)' INDEX EB/NO (DB) PROBABILITY PROBABILITY 


& BIT ERROR’ 
EBNODB = INIT 
DO 10 I=1,NUMPTS 


EBDVNO = 10. ODO**(EBNODB/10) 

SREDSI = DSQRT( 2*EBDVNO*K) 

Al = -SREDSI 

INTGL1 = DCADRE(F1,A1,B1,AERR,RERR,ERROR,IER) / SQR2PI 
INTGL2 = DCADRE(F2,A2,B2,AERR,RERR,ERROR,IER) / SQR8PI 
INTGL3 = DCADRE(F3,A1,B2,AERR,RERR,ERROR,IER) / SQR8PI 


HE CINTGES] Gi. 0;0D0) THEN 
PRBTCR = (N-1)*( INTGL2+INTGL3) + INTGL1 


ELSE 
PRBTCR = (N-1)*INTGL2 + (N-1)*INTGL3 + INTGL1 
libs 
PRBERR = ONE - PRBTCR 
DBDPBE = DERFC( DSQRT(2*K*EBDVNO/M) ) / 2 
UBAPBE = ( (M - 2)*DERFC( DSQRT(K*EBDVNO/2) ) 


& + DERFC( DSQRTCK*EBDVNO) ) ) /4 
WRITE( 8,101) I,EBNODB,PRBERR,DBDPBE , UBAPBE 
EBNODB = EBNODB + DLTADB 

10 CONTINUE 


100 FORMAT(2X,'NUMBER OF BITS: ',2X,11,6X,'NUMBER OF FREQUENCIES: ', 


& ei) 
101 FORMAT( 2X,13,3X,F10. 5,3X,E14. 5E2,3X,E£14. 5E2,3X,E14. 5E2) 
STOP 
END 


Ge POTS OTN Cas ee ae WUE Te COCOA TS ES REE ETCR TERT RR Res He Se fe sotesotedesede se setetedede sede detected ER EREE 


DOUBLE PRECISION FUNCTION F1(Y) 


DOUBLE PRECISION Y,SREDSI,SQRT2 
INTEGER N 


COMMON SREDSI,N 


C INITIALIZE CONSTANT 
SQRT2 = DSQRT(2. ODO) 
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C COMPUTE F1(Y,SREDSI ,N) 
Fl = DERF( (Y+SREDSI)/SQRT2 )**(N-1) * DEXP(-(Y**2)/2) 


RETURN 
END 


Sede de dede Te eae a se seve te ese Fea He Ie Fee TE Fede Tee Te Ie VOTE TET Ie TEE Tee VET FEE Te Tee TITER IC ete Ke KI KR KI 


DOUBLE PRECISION FUNCTION F2(Y) 
DOUBLE PRECISION Y,FACTR1,FACTR2,SREDSI 
DOUBLE PRECISION SQRT2 


INTEGER N 
COMMON SREDSI,N 
C INITIALIZE CONSTANTS 
SQRT2 = DSQRT(2. ODO) 
C COMPUTE F2(Y,SREDSI,N) 
FACTR1 = DERF( Y/SQRT2 )**(N-2) 
FACTR2 = DERF( (Y+SREDSI)/SQRT2 ) 
LEZ = FACTR1*FACTR2*DEXP( -( Y**2) /2) 
RETURN 
END 


keavdevedstetedcddedetesedescdestcded devevesesedee ack se sesetode sk de dea weve de sede desk a dese rete de Fete Teese Fee Ke ee eT Ne Ie sete 96 


DOUBLE PRECISION FUNCTION F3(Y) 
DOUBLE PRECISION Y,FACTR1,FACTR2,SREDSI 
DOUBLE PRECISION SQRT2 


INTEGER N 
COMMON SREDSI ,N 

C INITIALIZE CONSTANTS 
SQRT2 = DSQRT(2. 0D0) 

C COMPUTE F3(Y,SREDSI ,N) 
FACTR1 = DERF( (Y+SREDSI)/SQRT2 )**(N-2) 
FACTR2 = DERF( Y/SQRT2 ) 
F3 = FACTR1*FACTR2*DEXP(-CY+SREDSI)**2/2 ) 
RETURN 
END 
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APPENDIX  C. 


16-PSK RAYLEIGH FADING BER COMPUTER 


PROGRAM 


PEE SES SOE ORS ES IS GAS SE SS oe i ote Gogh shine isc risey seer ise vier Skis aridem casei ists aa cs che ioriaeiseiseistio cmos case tseiseapancaers 


PROGRAM TO COMPUTE THE PROBABILITY OF A BIT ERROR FOR 16-PSK USING 


ENVIROMENT. 


fy 


* 
ae THE DIRECT BIT DETECTION (DBD) RECEIVER IN THE RAYLEIGH FADING 
* 
* 


SESE RUSE ES ESTES CSS ST ad Saeed ae aint Fabia ar ior ena eto aoeiae toe ia el sripeia iar seigeisrises SeiarAoe iri stoeior deel seiseisel se bickselscisc toes 


PROGRAM DBDRAY 


IMPLICIT CHARACTER*1 (A-Z) 


DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
DOUBLE PRECISION 
INTEGER 


EXTERNAL 
COMMON 


Gob ERSET(S,0,-1) 


EBNODB, INIT, FINAL, DLTADB , PBERAY 
EBDVNO, AVESNR, ALPHA, ALSDEG , BETA 
P1BERY,V1SFTN(4) , ARGMNT 

P2BERY , V2SFTN(4) ,ARGS$1,ARGS$2 

P3BERY , INTGRL(8) , TERM$1, TERM$2, TERMS3 
F1,A,B,ERRABS , ERRREL, INTGL1,ERREST 
F2,F3,INTGL2, INTGL3 

PI ,ONE,SUM$1,SUM$2 

I, INDEX, NUMPTS , IRULE 


E243 
AVESNR, ALPHA, BETA, INDEX 


SUPPRESS IMSL10 ERROR MESSAGES. 


INITIALIZE CONSTANTS 


ea = 3.141592654D0 

ONE = 1.0D0 

ALSDEG = 14. 25D0 

ALPHA = ALSDEG * PI/180 

Pers = PI/4 - ALPHA 

V1ISFTN(1) = DCOS( ALPHA) 

V1SFIN(2) = DCOS( BETA) 

V1ISFTN(3) = DSIN(CALPHA) 

V1ISFTN(4) = DSIN( BETA) 

V2SFTN(1) = DCOS(ALPHA) - DSIN( ALPHA) 
V2SFTN(2) = DCOSCALPHA) + DSIN( ALPHA) 
V2SFTN(3) = DCOS(BETA) - DSIN( BETA) 
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aan 


10 


20 


30 


V2SFTN(4) = DCOS( BETA) + DSIN( BETA) 


NUMPTS = 101 

INIT = -10. 0D0 

FINAL = 20. 0D0 

DLTADB = (FINAL - INIT)/(NUMPTS-1) 
A = 0. 0D0 

B = PI/8 

ERRABS = 0. 0D0 

ERRREL = 0.000 

IRULE = 2 ; 


COMPUTE PROBABILITY OF BIT ERROR 
FOR VALUES OF EBNODB IN THE RANGE OF: 
INIT < EBNODB < FINAL 


WRITE(8,*)'EFFECT OF RAYLEIGH FADING ON 16-PSK, DBD RCVR' 
WRITE(8,100) ALSDEG 
WRITE(8,*)' ' 


WRITE(8,*)' DBD RCVR' 
WRITE(8,**)' BIT ERROR' 
WRITE(8,*)'INDEX EB/NO (DB) PROBABILITY’ 


EBNODB = INIT 


BEGIN MAIN LOOP. 
DO 70 I=1,NUMPTS 
EBDVNO 10. ODO**"( EBNODB/10) 
AVESNR 2 * EBDVNO 


CALCULATE PROBABILITY OF 1B ERROR 
SUMS1 = 0. ODO 
DO 10 INDEX=1,4 
ARGMNT = EBDVNO / (4 * EBDVNO * VISFTINCINDEX)**2 + ONE) 
SUM$1 = SUMS$1 + V1SFTN( INDEX) * DSQRT( ARGMNT) 
CONTINUE 


PIBERY = (2.0D0 - SUM$1)/4 


CALCULATE PROBABILITY OF 2B ERROR 
SUM$1 = 0. ODO 
DO 20 INDEX=1,4 
ARGS1 = EBDVNO / (EBDVNO * V2$FTN(INDEX)**2 + ONE/2) 
SUMS1 = SUM$1 + V2SFTN( INDEX) * DSQRT(ARG$1) 
CONTINUE 


SUM$2 = 0. ODO 
DO 30 INDEX=1,4 
ARGS$1 = DSQRT( V2SFTN(INDEX)**2 + ONE / (2*EBDVNO) ) 
ARGS2 = ARG$1 / V2SFTN( INDEX - (-1)**( INDEX) ) 
SUM$2 = SUM$2. + V2SFTN(INDEX) * DATAN(ARGS$2) / ARGS$1 
CONTINUE 


P2BERY = ONE/2 - SUMS1/4 + SUM$2/(2*PI) 


CALCULATE PROBABILITY OF 3B ERROR 
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DO 40 INDEX=1,8 
CALL DQDAG (F1,A,B,ERRABS ,ERRREL, IRULE, INTGL1,ERREST) 
TERM$1 = DSQRT(AVESNR)*INTGL1/4 
CALL DQDAG (F2,A,B,ERRABS ,ERRREL, IRULE, INTGL2, ERREST) 
TERM$2 = INTGL2 
CALL DQDAG (F3,A,B,ERRABS ,ERRREL, IRULE, INTGL3 ,ERREST) 
TERM$3 = DSQRT(AVESNR)*INTGL3 
INTGRL( INDEX) = TERM$1 - (TERM$2 + TERM$3)/(4*PI) 

40 CONTINUE 


SUM$1 = 0. ODO 
DO 50 INDEX=5,8 
SUM$1 = SUM$1 + INTGRL( INDEX) 
50 CONTINUE 


SUM$2 = 0. 0D0 
DO 60 INDEX=1,4 
SUM$2 = SUM$2 + INTGRL( INDEX) 
60 CONTINUE 


P3BERY = ONE/2 + SUMS1 - SUMS2 
PBERAY = PIBERY/2 + P2BERY/4 + P3BERY/4 
WRITE( 8,101) I,EBNODB, PBERAY 
EBNODB = EBNODB + DLTADB 
70 CONTINUE 
C END MAIN LOOP 


100 PORMAT( 2X, ‘ALPHA =" ,1X,F5.2,1X,' DEGREES. °) 
Moe PORMAT( 2X,13,3X,F10.5,3X,614. 5E2,3X,E14. 5E2,3X,E14. 5E2) 


STOP 
END 


Fee Te Tete eee Tee Gee Te STE TE TS ESTE TS TS TS ESTES ESTEE TEES TET TESTS ETS TEES TECTED TETEEE TELE TE TEEPE EEE EE, 
DOUBLE PRECISION FUNCTION F1(GAMMA) 
IMPLICIT CHARACTER*1 (A-Z) 
DOUBLE PRECISION GAMMA,AVESNR,ALPHA, BETA 
DOUBLE PRECISION TERM$1,FINDEX,GINDEX 
INTEGER INDEX 


COMMON AVESNR , ALPHA, BETA , INDEX 


E COMPUTE F1( INDEX ,GAMMA, AVESNR, ALPHA, BETA) 
GO TO (10,20,30,40,50,60,70,80) INDEX 


10 FINDEX = DSIN(GAMMA - ALPHA) 
GINDEX = DCOS(GAMMA - ALPHA) 
GO TO 100 
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20 FINDEX = DCOS(GAMMA - ALPHA) 
GINDEX = DSIN(GAMMA - ALPHA) 
GO TO 100 

30 FINDEX = DCOS(GAMMA + ALPHA) 
GINDEX = DSIN(CGAMMA + ALPHA) 
GO TO 100 

40 FINDEX = DSIN(GAMMA + ALPHA) 


GINDEX = DCOS(GAMMA + ALPHA) 
GO TO 100 


50 FINDEX = DSIN(GAMMA - BETA) 
GINDEX = DCOS(GAMMA - BETA) 


GO TO 100 

60 FINDEX = DCOS(GAMMA - BETA) 
GINDEX = DSIN(GAMMA - BETA) 
GO TO 100 

70 FINDEX = DCOS(GAMMA + BETA) 
GINDEX = DSIN(GAMMA + BETA) 
GO TO 100 

80 FINDEX = DSIN(GAMMA + BETA) 
GINDEX = DCOS( GAMMA + BETA) 

100 TERM$1 = AVESNR * FINDEX**2 + 1.0D0/2 
Veal = GINDEX / TERMS1 / DSQRT(TERMS1) 
RETURN 
END 


Petede dete de teds Tete ds ce Tete Te Tee Te TOTS TS TEES TTT EE EEE TET ETE EEL EE TEER EEE EERE REREEERRRBREE 
DOUBLE PRECISION FUNCTION F2(GAMMA) 
IMPLICIT CHARACTER*1 (A-Z) 
DOUBLE PRECISION GAMMA,AVESNR ,ALPHA, BETA 
DOUBLE PRECISION TERM$1,TERMS2,FINDEX,GINDEX 
INTEGER INDEX 


COMMON AVESNR , ALPHA, BETA, INDEX 


fe COMPUTE F2( INDEX,GAMMA, AVESNR , ALPHA, BETA) 
GO TO (10,20,30,40,50,60,70,80) INDEX 
10 FINDEX = DSIN(GAMMA - ALPHA) 
GINDEX = DCOS(GAMMA - ALPHA) 
GO TO 100 


20 FINDEX = DCOS(GAMMA - ALPHA) 
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GINDEX = DSIN(GAMMA - ALPHA) 
GO TO 100 


30 FINDEX = DCOS(GAMMA + ALPHA) 


GINDEX = DSINCGAMMA + ALPHA) 
GO TO 100 

40 FINDEX = DSIN(GAMMA + ALPHA) 
GINDEX = DCOS(GAMMA + ALPHA) 
GO TO 100 

50 FINDEX = DSIN(GAMMA - BETA) 
GINDEX = DCOS(GAMMA - BETA) 
GO TO 100 

60 FINDEX = DCOS(GAMMA - BETA) 
GINDEX = DSIN(GAMMA - BETA) 
GO TO 100 

70 FINDEX = DCOS( GAMMA + BETA) 
GINDEX = DSIN(GAMMA + BETA) 
GO TO 100 

80 FINDEX = DSIN(GAMMA + BETA) 
GINDEX = DCOS( GAMMA + BETA) 

HOOT TERMS I = AVESNR * FINDEX**2 + 1: 0D0/2 
TERMS2 = AVESNR * GINDEX**2 
BZ =senis2 /{ CIERMSI-TERMS2 + TERHS1** 2) 
RETURN 
END 


DOUBLE PRECISION FUNCTION F3(GAMMA) 

IMPLICIT CHARACTER*1 (A-Z) 

DOUBLE PRECISION GAMMA,AVESNR,ALPHA,BETA 

DOUBLE PRECISION TERM$1,TERMS2,FACTR1,FACTR2 ,F INDEX ,GINDEX 
INTEGER INDEX 


COMMON AVESNR , ALPHA , BETA , INDEX 


C COMPUTE F3( INDEX ,GAMMA,AVESNR, ALPHA , BETA) 


GO TO (10,20,30,40,50,60,70,80) INDEX 


10 FINDEX = DSINCGAMMA - ALPHA) 
GINDEX = DCOS(GAMMA - ALPHA) 
GO TO 100 


20 FINDEX = DCOS(GAMMA - ALPHA) 
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30 


40 


50 


60 


70 


80 


100 


GINDEX = DSIN( GAMMA 
GO TO 100 


ALPHA ) 


+ 


1. 0D0/2 


FINDEX = DCOS(GAMMA + ALPHA) 
GINDEX = DSIN(GAMMA + ALPHA) 
GOTO 100 

FINDEX = DSIN(GAMMA + ALPHA) 
GINDEX = DCOS(GAMMA + ALPHA) 
GO TO 100 

FINDEX = DSIN(GAMMA - BETA) 
GINDEX = DCOS(GAMMA - BETA) 
GO TO 100 

FINDEX = DCOS(GAMMA - BETA) 
GINDEX = DSIN(GAMMA - BETA) 
GO TO 100 

FINDEX = DCOS(GAMMA + BETA) 
GINDEX = DSIN(GAMMA + BETA) 
GO TO 100 

FINDEX = DSIN(GAMMA + BETA) 
GINDEX = DCOS(GAMMA + BETA) 
TERMS1 = AVESNR * FINDEX**2 
TERM$2 = AVESNR * GINDEX**2 
FACTR1 = GINDEX / TERM$1 / DSQRT(TERMS1) 
FACTR2 = DATAN( DSQRT(TERMS2/TERMS1) ) 
F3 = FACTR1 * FACTR2 
RETURN 

END 
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APPENDIX  D. MFSK/QPSK RAYLEIGH FADING BER COMPUTER 


PROGRAM 
PES LS LE SAA Safir c fo ala cose sos tor fede beri tse dee toleari el sei set srioid sec teige lack acicel ashi seasehseaschociscks 
oS 
oye 
ve PROGRAM TO COMPUTE THE BIT ERROR RATES FOR THE BI-ORTHOGONAL 
* SIGNAL SET (BOSS) FOR THE OPTIMUM RCVR WHEN SUBJECTED TO RAYLEIGH 
* FADING. 


GETCTTCSCG eae WISTS Te PETE TC TC TOTS TCCR TCT IC EC CS TC TE TCTCTE FET CTOTO TOUT CTE TCT TE TON TCT, TOA TCTSTOTOTO TOTTI IC TOP RTT 


PROGRAM PBTERA 


IMPLICIT CHARACTER*1 (A-Z) 

DOUBLE PRECISION EBNODB,INIT,FINAL,DLTADB 

DOUBLE PRECISION EBDVNO,SOSQRD,S 

DOUBLE PRECISION PBCRAY,PBERAY 

DOUBLE PRECISION DCADR2,INTGND,LIL,UIL, AERR,RERR, ERROR 
DOUBLE PRECISION ONE 


INTEGER I,K,N,M,F, IER,NUMPTS 
INTEGER LEVEL, LEVOLD 
EXTERNAL INTGND 
COMMON SOSQRD,S,N 

€ INITIALIZE CONSTANTS 
K = 5 
M — Dek K 
N = Pao K- |) 

F = 2**(K-2) 

LiL = 00. ODO 

ONE = 01. 0D0 

NUMPTS = 101 

INIT = -10.0D0 

FINAL = 20. 0D0 

DLTADB = (FINAL - INIT)/(NUMPTS-1) 
AERR = 0.0 

RERR = .0001 

C SUPPRESS ALL ERROR MESSAGES. 
LEVEL = 0 
CALL UERSET( LEVEL, LEVOLD) 

C CALL XUFLOW(0) 
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COMPUTE PROBABILITY OF BIT ERROR 
FOR VALUES OF EB/NO IN THE RANGE OF: 
INIT < EBNODB < FINAL 


aan 


WRITE(8,*)'EFFECT OF RAYLEIGH FADING ON BOSss' 
WRITE(8,*)' ' 
WRITE(8,100) K,F 


WRITE(8,*)' ' 

WRITE(8 ,*)' OPTIMUM RCVR' 
WRITE(8,**)' BIT ERROR ' 
WRITE(8,*)' INDEX EB/NO (DB) PROBABILITY ' 


EBNODB = INIT 
DO 10 I=1,NUMPTS 
EBDVNO = 10. OD0**(EBNODB/10) 


SOSQRD = EBDVNO * K 
UIL = MIN(18.5D0,10*DSQRT(SOSQRD) ) 
PBERAY = DCADR2( INTGND,LIL,UIL,AERR,RERR,ERROR,IER) / SOSQRD 


WRITE(8,101) I,EBNODB,PBERAY 
EBNODB = EBNODB + DLTADB 
10 CONTINUE 


100 FORMAT(2X,'NUMBER OF BITS: ',2X,11,6X,'NUMBER OF FREQUENCIES: ', 
& 2X, aii) 
101 FORMAT( 2X,13,3X,F10.5,3X,E14.5E2, 3x, E14, Seon ees 


STOP 
END 


aededededevedodedede ceded dedevodevede decked ke sesevedeveved seve sevevede se ses sek Ketetevedesek Kester He Ketek se ve tek Fete sete He Hk eH Fe 
DOUBLE PRECISION FUNCTION INTGND( SREDSI) 


IMPLICIT CHARACTER*1 (A-Z) 

DOUBLE PRECISION EBDVNO,SREDSI,SOSQRD,S,THRHLD 
DOUBLE PRECISION PRBTCR,PRBTER,ORTERM,APTERM 
DOUBLE PRECISION INTGL1,INTGL2,INTGL3 

DOUBLE PRECISION DCADRE,F1,A1,B1,AERR,RERR,ERROR 


DOUBLE PRECISION F2,A2,B2,F3 
DOUBLE PRECISION PIE,SQR2PI,SOQR8PI,ONE 
INTEGER N, IER 
EXTERNAL Bign2 ous 
COMMON SOSQRD,S,N 

c INITIALIZE CONSTANTS 
THRHLD = 6. 2D0 
PIE = 3.141592654D0 
ONE = 1.00D0 
SQR2PI = DSQRT(2*PIE) 

SQR8PI = DSQRT(8*PIE) 
AERR = 0.0 
RERR = .0000001 


OZ 


A2 = 00. ODO 
Bl = 20. 0D0 
B2 = 20. 0D0 
C COMPUTE PROBABILITY OF BIT ERROR AS A 
é FUNCTION OF SREDSI AND N 
IF(SREDSI.GT. THRHLD) THEN 
G USE UNION BOUND APPROXIMATION 
ORTERM = (2*N - 2)*DERFC( SREDSI/2 ) 
APTERM = DERFC( SREDSI / DSQRT(2.0D0) ) 
PRBTER = ( ORTERM + APTERM ) / 4 
ELSE 
C USE EXACT CALCULATION 
S = SREDSI 
Al = -SREDSI 
INTGL1 = DCADRE(F1,A1,B1,AERR,RERR,ERROR,IER) / SQR2PI 
INTGL2 = DCADRE(F2,A2,B2,AERR,RERR,ERROR,IER) / SQR8PI 
INTGL3 = DCADRE(F3,A1,B2,AERR,RERR,ERROR,IER) / SQR8PI 


MECINGGL3.60. 0. 0D0) THEN 


PRBTCR = INTGL1 + (N-1)*( INTGL2+INTGL3) 
ELSE 
PRBTCR = INTGL1 + (N-1)*INTGL2 + (N-1)*INTGL3 
ENDIF 
PRBTER = ONE - PRBTCR 
ENDIF 
C COMPUTE INTGND(SREDSI,SOSQRD,N). 


INTGND = PRBTER * SREDSI * DEXP( -(SREDSI**2) / (2 * SOSQRD) ) 


RETURN 
END 


Sodededsdods ton ants JeJanleds dace n'ewsdewdes Sodas JentenlededacladedededodensSedodadodetervledotaderisdatadededevlededodode dade dedeSodedawdeclectedadadau'sJaclauis 
POOC ES SCC GC BN ee aR aN AN EC EN EN EN EN EN EVES EDEN ES ELEN ECTS ECU EC ESE EC EN EC ECTS ENE EC OCIS ES EC EC EN TSEC ES EC EC ER ECE TC ES EN ES ES EC ES EVENTS ES EC ES EC ENED ENED OD 


DOUBLE PRECISION FUNCTION Fi(Y) 


DOUBLE PRECISION Y,SOSQRD,S,SQRT2 


INTEGER - N 
COMMON SOSQRD ,S,N 

C INITIALIZE CONSTANT 
SQRT2 = DSQRT(2. ODO) 

C COMPUTE FIUCY,S,N) < 
Fl = DERF( (Y+S)/SQRT2 )**(N-1) * DEXP(-(Y**2)/2) 
RETURN 
END 
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Feredesedkse 


Seitcdetededs 


dededededede de deve ded dete dete de de deve Keve dete ae Fe deve He de det He dete sede de ete eee de Fee eK RK Kae eK 


DOUBLE PRECISION FUNCTION F2(Y) 
DOUBLE PRECISION Y,FACTR1,FACTR2,SOSQRD,S 
DOUBLE PRECISION SQRT2 


INTEGER N 
COMMON SOSQRD,S,N 

INITIALIZE CONSTANTS 
SQRT2 = DSQRT(2. 0D0) 


COMPUTE F2(Y,S,N) 


FACTR1 = DERF( Y/SQRT2 )**(N-2) 

FACTR2 = DERF( (Y+S)/SQRT2 ) 

F2 = FACTR1*FACTR2*DEXP( -(Y**2)/2) 
RETURN 

END 


detedetedetetedede dete dete teds Tete ve Te de sek Fete dese se te Te Fe Fe TT TTF TEN TITER IRE RICK TIER 
DOUBLE PRECISION FUNCTION F3(Y) 

DOUBLE PRECISION Y,FACTR1,FACTR2,SOSQRD,S 

DOUBLE PRECISION SQRT2 


INTEGER N 
COMMON SOSQRD,S,N 

INITIALIZE CONSTANTS 
SQRT2 = DSQRT(2. ODO) 


COMPUTE ESCY ¢5...N) 


FACTR1 = DERF( (Y+S)/SQRT2 )**(N-2) 
FACTR2 = DERF( Y/SQRT2 ) 

F3 = FACTR1*FACTR2*DEXP( -(Y+S)**2/2 ) 
RETURN 

END 
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